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Abstract. The debate on the oxygen abundances of metal-poor stars has its origin in contradictory results ob¬ 
tained using different abundance indicators. To achieve a better understanding of the problem we have acquired 
high quality spectra with the Ultraviolet and Visual Echelle Spectrograph at VLT, with a signal-to-noise of the 
order of 100 in the near ultraviolet and 500 in the optical and near infrared wavelength range. Three different 
oxygen abundance indicators, OH ultraviolet lines around 310.0 nm, the [O i] line at 630.03 nm and the OI lines 
at 777.1-5 nm were observed in the spectra of 13 metal-poor subgiants with —3.0 < [Fe/H] < —1.5. Oxygen abun¬ 
dances were obtained from the analysis of these indicators which was carried out assuming local thermodynamic 
equilibrium and plane-parallel model atmospheres. Abundances derived from OI were corrected for departures 
from local thermodynamic equilibrium. Stellar parameters were computed using T e g-vs-color calibrations based 
on the infrared flux method and Balmer line profiles, Hipparcos parallaxes and Fell lines. [O/Fe] values derived 
from the forbidden line at 630.03 nm are consistent with an oxygen/iron ratio that varies linearly with [Fe/H] as 
[O/Fe] = —0.09(±0.08)[Fe/H] + 0.36(±0.15). Values based on the Oi triplet are on average 0.19 ± 0.22dex(s.d.) 
higher than the values based on the forbidden line while the agreement between OH ultraviolet lines and the 
forbidden line is much better with a mean difference of the order of —0.09 ± 0.25 dex(s.d.). In general, our results 
follow the same trend as previously published results with the exception of the ones based on OH ultraviolet lines. 

In that case our results lie below the values which gave rise to the oxygen abundance debate for metal-poor stars. 

Key words, line: formation - stars: abundances - stars: atmospheres - stars: population II - Galaxy: evolution 


1. Introduction 

Oxygen plays a major role in astrophysics through its 
high cosmic abundance. It is predominantly produced and 
ejected back into the interstellar medium in connection 
with core-collapse supernovae (SNe II). By comparing oxy¬ 
gen abundances in stars of different ages with those of iron, 
which is produced both in SNe II and in thermo-nuclear 
supernovae of white dwarfs (SNe la), one can probe the 
star formation rate history and initial mass funct ion of a 
galaxy as well as the phy sics of supernovae (e.g. iTinslevI 
fl980nWheeler et all198911 . Due to the longer time-scale of 
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Fe production in SNe la , one can expect an oxygen over¬ 
abundance relative to iron at low metallicity: [O/Fe] > 0 1 . 

While such an oxygen over-abund ance in metal-poo r 
stars has long been known to exist ( Conti et HD Em 
the exact amount is still contested with no unanimous 
agreement apparently in sight. Broadly speaking the de¬ 
rived [O/Fe] values depend on which oxygen diagnos¬ 
tic is employed in the analysis and what type of stars 
is observed. There are four different types of spectral 
transitions which have been utilised for this purpose: 
the forbidden [Oi] 630.0 and 636.3nm lines, high exci¬ 
tation Oi lines, in particular the Oi 777.1-5nm triplet, 
OH A-X electronic lines in the ultraviolet (UV) and OH 
vibration-rotation lines in the infrared (IR). The Oi and 
OH lines have mainly been used in metal-poor turn-off 
stars (T e ff^ 6000 K) while the O abundances in halo gi¬ 
ants (T e g < 5000K) are typically based on the [Oi] lines. 


1 The abundance ratios are defined by the customary 

[X/Fe] = log (Nx/We)* - log (N x /N Fe ) Q . 












2 


A.E. Garcia Perez et al.: Oxygen abundances in metal-poor subgiants 


The [Oi] lines suggest a quasi-plateau at [O/Fel ~ +0-5 
for [Fe/H] < —1 (e.g. lBarbuvlll988l: ISneden et al.lll99lt 
INissen et aD^^). This has traditionally been advocated 
as the correct metallicity trend, in particular since the first 
studies based on the OH l ines in the UV gave con sistent 
results ( Bessell et al.ll984ll99itlNissen et all 1994. More 
recently, the OH lines in the IR have also been studied in 
a few stars with a similar result, at l east when restricting 
to weak, higher excitation OH lines llBalachandran et all 
2001 : iMelendez et al . 2 Qq 3). Systematically higher [O/Fe] 
value s are normally derive d when using the OI triplet 
fe.g. lAbia fc Rebold Il989t llsraelian et a,lJ Il998l l 200 ll 
Bocsgaard et al.lll999t ICarretta et alJl 200 (I INissen et al.l 
20021 Fulbright fc .Iohnsonll2003(l . sometimes interpreted 
as a steady linear increase in [O/Fe] towards lower [Fe/H]. 
The OI triplet results have often been blamed on depar¬ 
tures from local thermodynamic equilibrium (LTE), ef¬ 
fects of inhomogeneities and/or an erroneous T e g-scale in 
order to reconcile them with the [O i] results. The issue of 
oxygen abundances in halo stars received much renewe d 
attention fo llowi ng the surprising results of llsraelian et all 
lll998H200lh and lBoesgaard et a,l.l ill 999ll who found a near 
linear trend in [O/Fe] with a significant slope of about 
—0.4 towards lower metallicity from OH lines in the UV 
in turn-off stars and a couple of subgiants. 

All of the above-mentioned oxygen diagnostics have 
their advantages and disadvantages. The O I tri p let is sus¬ 
ceptible to non-LTE effects (NLTE) llKiselrna.nll200l[ and 
references therein) and is, due to the high excitation po¬ 
tential, quite sensitive to the adopted T e g. The [Oi] line 
is immune to departures from LTE but is essentially un¬ 
detectable in turn-off stars with [Fe/H] ^ —2 due to its 
weakness. The OH lines are very sensitive to the temper¬ 
ature structure in the stellar atmosphere. Given the much 
lower temperatures encountered in the optically thin lay¬ 
ers in 3D hydrodynamical model atmospheres compared 
with classical ID hyd rostatic models for metal-poor stars 
iAsplund et alJll999 : h large downward abundance correc¬ 
tions have been flagged for the OH lines in turn-off stars. 
The [Oi] lines are also affected by such temperature in¬ 
homogeneities but to a lesser degree than the OH lines 
iNissen et al.1 Eo02 h while the NLTE line formation of 
the Oi lines can be suspected to be depen d ent on the 
3D atmospheric structure as well iKiselman &; Nor dlundl 

1995; Asplund fe Garcia PereJ 200 lHAsplund et al.l 2QQ4|) . 

Additional complications and confusion stem from the dif¬ 
ferences and uncertainties in the adopted T e g and log < 7 , 
the use of Fe 1 or Fe 11 lines to derive [Fe/Hl, poss ib le miss¬ 
ing continuous opacit ies in the UV l|Balachandran fc Belli 
ll 998t iBell et aT1l200lh and which solar oxygen abundance 
the stellar values are referenced to. Given the recent large 
dow nward revision of the solar O abundance to log eo ~ 
8.7 l|Allende Prieto et al .112(10 ll I Asplund et al .11200411 . this 
difference of abou t 0.2 dex compared with t he previously 
advocated value ( Anders & GrevesseiIl989 1 will directly 
translate to a corresponding increase in [O/Fe]. 

Considering all these remaining uncertainties it is per¬ 
haps not surprising that no consensus has as yet been 


achieved in terms of the behaviour of [O/Fe] with [Fe/H]. 
It is noteworthy that significant differences between the 
different oxygen indicators only reveal themselves at low 
metallicities ([Fe/H] <) —2), which is observationally very 
challenging due to the weakness of some of the spec¬ 
tral features. The fact that the different oxygen diagnos¬ 
tics are normally used in different types of stars consti¬ 
tutes a major problem in this regard. The only metal- 
poor stars for which it is possible to employ molecular 
as well as forbidden and permitted atomic lines are sub¬ 
giants but also then exceptionally high-quality spectra are 
required. Until now, only two halo subgiants (HD 140283 
and BD+23°3130, both with [Fe/H] ~ — 2.4) have been 
analysed wi t h [Oil, Oi and OH UV lines dlsraelianet_a]J 
199R 200ll: iFulbright fc Krafd ll99H ICavrel et all l200ll: 


ign_ _ _ _ _ _ — 

Balachandran et alJl200lt ~ fNissenet alJl2002|h with some¬ 

what ambiguous results. Our best hope to finally resolve 
this important outstanding problem is offered by halo sub¬ 
giants. In the present article we present such a study based 
on very high-quality UV and optical spectra for a sample 
of 13 subgiants with —3.0 < [Fe/H] < —1.5. 


2. Observations 

We identified a sample of halo subgiants based on 
Hipparcos parallaxes and/or Stromgren photometry with 
metallicity —3.0 < [Fe/H] < —1.5. From this sample, 
high-quality sp ectr a were ac quired for 13 subgiants us¬ 
ing the UVES llPekker et al l 20001 at ESO’s Very Large 
Telescope (VLT) 8 m Kueyen telescope. The observations 
were carried out in service mode between October 2000 
and March 2001. The use of the dichroic capability of 
UVES enabled simultaneous observations of the UV and 
optical region to cover the three most important oxygen 
abundance diagnostics: OH A-X lines around 310.0 nrn , 
the [Oi] 630.03nm line and the Oi 777.1-5nm lines. Due 
to the relatively low effective temperature of the targets, 
the observing times were largely determined by the re¬ 
quirement to obtain reasonable signal-to-noise ( S/N ) in 
the UV (resulting in S/N ~ 100-150 at 310.0 nm ) . This 
enabled the total exposure time to be split between two 
optical settings to cover both the forbidden and permitted 
0 1 lines and still achieve very high S/N (> 500 in the cen¬ 
tre of an echelle order at both 630 and 777 nm, although 
as will be described further below the 01 triplet falls close 
to the edge of two orders in most cases resulting in typical 
values of S/N ~ 250 at the triplet). We thus used two dif¬ 
ferent dichroic settings for the red arm of the spectrograph 
(580 and 860 nm) together with the 346 nm setting for the 
blue arm. In total, the spectra covered 302.4-388.4nm and 
476.4-1061.Onm, with near complete coverage in these re¬ 
gions. In all cases, at least three different exposures were 
obtained in each setting to allow efficient removal of cos¬ 
mic rays. 

For the blue arm, a projected slit-width of 1.0" was 
used while for the red arm the width was 0.7". The re¬ 
sulting nominal spectral resolving powers associated with 
these slit-widths were 45 000 and 55 000, respectively, with 
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4 pixels per spectral resolution element. In order to gain 
S/N in the UV and because of oversampling, a 2x2 bin¬ 
ning was used when reading out the EEV CCD without 
loosing any spectral resolution. 

As there are several telluric lines due to in particu¬ 
lar H 2 0 and 0 2 around 630.03 nrn, the [Ol] line could 
be blended by telluric features depending on the radial 
velocity of the star and the particular date of the ob¬ 
servations. We tried to avoid such problems by placing 
time constraints on when the observations should be car¬ 
ried out in service mode. For this purpose, we obtained 
radial velocity information of our targets from SIMBAD 
and computed the location of the stellar line in relation 
to the telluric lines for the relevant observing period. The 
procedure worked very well with one exception in which 
the radial velocity reported by SIMBAD turned out to 
be different from the actual value (given in Tabled by 
about 20kms _1 . As a result, our spectra of this star 
(CD —24° 1782) shows the [O i] 630.03 nrn significantly dis¬ 
torted by telluric absorption, unfortunately rendering it 
useless for an abundance analysis. 

The procedure followed in the data reduction was the 
classical one: subtraction of bias, definition of the echelle 
orders, division by the normalised flat-field, extraction of 
the orders, wavelength calibration and continuum normal¬ 
ization, all carried out within iraf. The continuum was 
fitted using cubic spline functions. Due to the metal-poor 
nature of our targets, it was always relatively straight¬ 
forward to identify wavelength regions without significant 
line contributions in the optical to make the continuum 
well defined. The UV region is more crowded even in 
metal-poor stars but it was still possible to find suffi¬ 
ciently many apparently clean continuum regions in the 
neighborhood of the OH lines. For the UV spectral re¬ 
gion, the S/N was modest and hence the UVES pipeline 
reduction worked well, making a special data reduction 
unwarranted. 

As mentioned above, care was taken when preparing 
for the observations in order to minimise the risk for blend¬ 
ing of the [O i] line with telluric lines. Additional precau¬ 
tions in this regard were taken by obtaining spectra of 
rapidly rotating B stars during the same nights and at 
about the same airmass as the programme stars. In a few 
cases when the stellar [Ol] line was partially blended by 
such telluric absorption, we relied on the iraf task tel¬ 
luric, which worked well. The task optimises the fit be¬ 
tween the telluric lines of the programme star and the B 
star by allowing a relative scaling in airmass and a shift in 
wavelength. The stars in our sample with partially blended 
[Oi] lines are HD 4306, HD 26169, HD 45282 , HD 126587, 
HD 128279 and HD 218857. 

Observed stellar radial velocities for the stars were de¬ 
rived from the Doppler shifts of three Cal (612.22, 616.22, 
643.91 nm) and one Fei (623.07nm) lines. The resulting 
heliocentric velocities are given in Tabled 


3. Model atmospheres and stellar parameters 

Before proceeding to determining element abundances 
for our sample, the stellar fundamental parameters must 
be estimated. The effective temperature (Teff), surface 
gravity (logg), metallicity ([Fe/H]) and microturbulence 
(Cturb) are the parameters needed to specify standard ID 
model atmospheres. Such model atmospheres are used to 
calculate theoretical stellar spectra which can be com¬ 
pared with observed spectra directly, or indirectly by com¬ 
paring observed and computed equivalent widths of spec¬ 
tral lines. 


3.1. Model atmospheres 


The transport of radiation in the stellar atmospheres 
needs to be modelled in order to enable comparison with 
the observed spectra. For the purpose, we have calculated 
theoretical ID LTE model atmospheres in hydro static 
equil i brium using the MA RCS program I Clnsta.fsson et ID 
ll 975tfAsnlimd et al. 1179 971 for the stellar parameters asso¬ 
ciated with the observed stellar sample. These models are 
based on up-to-date continuous opacities and include the 
effects of line-blanketing through opacity sampling. In the 
program, the convective energy flux is estimated using the 
mixing-length theory, assuming a mixing-length parame¬ 
ter a = 1.5. The models have all been computed using 
a microturbulence £turb = l.Okms -1 . The adopted solar 
element abundances f or the construction of the model at¬ 
mospheres stem from Idrevesse V Sanvall ( 1998 1 with the 
exception of O which is taken from the ID estimate of 
iNissen et alJ ([200211 : logeo= 8.74. The chemical compo¬ 
sitions of the models corresponding to our stellar sample 
are then scaled by the stellar metallicity, assuming that all 
the a-elements (with carbon not included) were enhanced 
by 0.4 dex. 

The solar oxygen and iron abundances used to estimate 
[O/Fe] and [Fe/H] were derived using a solar MARCS model 
atmosphere with the stellar parameters T e ff = 5780 K, 
logg = 4.44 [cgs] and [Fe/H] = 0.00 and adopting a micro¬ 
turbulence £turb = 1.15kms _1 for the spectral line forma¬ 
tion calculations. The basic reason for our choice of a theo¬ 
retical MARCS model to represent the solar atmosphere in¬ 
stead of a semiempirical one, such as the Holweger-Mueller 
model, was the ambition to make the analysis differential 
as far as possible, where systematic errors in models used 
for the Sun and the programme stars hopefully could can¬ 
cel to a considerable extent. 


3.2. Reddening 

The majority of stars in our sample are more distant than 
200 pc according to their Hipparcos parallaxes. For these 
distances interstellar extinction could become important 
and the observed stellar magnitudes and stellar colours 
have to be corrected of reddening. Fortunately, most of the 
stars are found at high galactic latitudes where reddening 
in general is limited. Two different methods have been con- 
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Table 1 . List of observed stars including coordinates, uvbyfl photometry, (V — K ) colour, Hipparcos parallaxes and 
the measured radial velocity from our observations 


Star 

a(2000) 

5(2000) 

V 

(b-y) 

mi 

Cl 

0 

(V - K) 

TV 

o-(tt) 

Vr 


[hr:min:s] 



1 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mas] 

[mas] 

[kins -1 ] 

HD 4306 

00:45:27.2 

-09 

: 32 : 

39.8 

9.035 

0.518 

0.052 

0.348 

2.529 

2.212 

4.81 

1.40 

-69.8 

HD 26169 

04:00:52.4 

-75 

: 36 : 

11.5 

8.797 

0.523 

0.051 

0.328 

2.521 

2.163 

2.83 

0.79 

-35.6 

HD 27928 

04:22:55.1 

-37 

: 15 : 

49.2 

9.554 

0.506 

0.067 

0.295 

2.523 

2.106 

2.36 

1.16 

43.3 

HD 45282 

06:26:40.8 

03 

: 25 : 

29.8 

8.028 

0.451 

0.108 

0.277 

2.544 

1.939 

7.34 

0.96 

305.7 

HD 108317 

12:26:36.8 

05 

: 18 : 

09.0 

8.036 

0.447 

0.055 

0.291 

2.548 

1.883 

4.53 

1.06 

6.6 

HD 126587 

14:27:00.4 

-22 

: 14 : 

39.0 

9.119 

0.602 

0.042 

0.460 


2.451 

1.40 

1.44 

149.6 

HD 128279 

14:36:48.5 

-29 

: 06 : 

46.6 

8.039 

0.467 

0.055 

0.262 

2.545 

1.974 

5.96 

1.32 

-76.0 

HD 200654 

21:06:34.7 

-49 

: 57 : 

50.3 

9.097 

0.460 

0.027 

0.271 

2.534 

1.943 

3.20 

1.25 

-46.1 

HD 218857 

23:11:24.6 

-16 

: 15 : 

04.0 

8.967 

0.506 

0.076 

0.371 


2.094 

3.51 

1.42 

-170.9 

HD 274939 

05:33:18.0 

-47 

: 56 : 

13.8 

9.435 

0.494 

0.128 

0.300 

2.526 

2.136 

2.88 

1.06 

172.0 

BD—01° 2582 

11:53:37.3 

-02 

: 00 : 

36.7 

9.595 

0.493 

0.097 

0.298 


1.997 

2.98 

1.35 

0.8 

CD—24° 1782 

03:38:41.5 

-24 

: 02 : 

50.3 

9.934 

0.466 

0.040 

0.287 

2.526 

1.975 

4.42 

1.75 

101.0 

CD—30° 0298 

00:58:43.9 

-30 

: 05 : 

57.7 

10.80 

0.483 

0.025 

0.296 


1.973 

1.06 

2.23 

27.2 


sidered for estimating redd ening: the maps of infrared dust 
emission bv ISchlegel et alj il998l) and the models of visual 
interstellar extinction bv lHakkila et al ] ill 997h . In addi¬ 
tion, the calibration of th e interstellar doub l et Na i lines 
at 589.0 and 589.6 nm bv lMunari Zwitt.eil ill 997tl have 
been tak en into acc ount in some cases for comparison pur¬ 
poses. ISchlegel et alJ lll998ll used the COBE/DIRBE and 
IRAS/ISSA infrared maps of dust emission over the entire 
sky to estimate the dust column density. To proceed from 
the dust column density to reddening they used colours of 
elliptical galaxies. An important aspect of t hese maps is 
the po ssible presence of filamentary details. lHakkila et all 
(19971 calculated reddenings based on published results of 
large-scale visual interstellar extinction. While the redden¬ 
ing estimates based on the model by Hakkila et al. depend 
on the assumed stellar distances, the estimates from the 
IR dust emission maps of Schlegel et al. give the reddening 
that an object will have if it would lie outside the dusty 
Galaxy. 

The estimated reddening values E( b-v) f° r both the 
Schlegel et al. and Hakkila et al. methods are given in 
Tabled together with the galactic longitude and latitude. 
The values based on Schlegel et al. are almost always 
higher than those derived from Hakkila et al. In Fig.m 
the differences between the two methods are displayed 
against stellar distance and galactic latitude. There may 
be an indication that the greatest differences occur for low 
galactic latitudes. It is difficult to assess which values are 
more accurate. The uncertainties in the reddening given 
by Schlegel et al. are claimed to be of the order of 16%, 
Hakkila et al. estimate the errors to be about 0.07 mag. 
Indeed, the differences we have found could be interpreted 
as an indication of the uncertainty in the estimated red¬ 
dening values. However, we believe that the Hakkila-based 
values may give a better description of the interstellar dust 
distribution in the solar neighborhood. We also note that 
the IR dust maps used by Schlegel et al. may be contami¬ 
nated by infrared point sources that could not be properly 
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Fig. 1. Compari son of reddening e stimates based on the 
two m ethods bv lHakkila et al J dl997frl and ISchlegel et all 
(Il998h . The differences are plotted versus the stellar dis¬ 
tance in parsecs (d) in the upper panel and versus the 
Galactic latitude of the star in degrees (b) in the bottom 
panel. Note that HD 45282 is not shown in the figure as it 
would fall well outside the plot limits (AI2( B _ V ) = 0.78). 


eliminated. One such example is HD 45282 for which the 
Schlegel et al. value is If(B_v) = 0.82. HD45282 is a rela¬ 
tively nearby star (d ~ 140 pc) so the IR dust emission 
maps most probably overestimate the reddening, although 
we note that the star is located at low latitudes. As will 
be seen below, with £ , ( B _v)=0.82 the effective tempera- 
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Table 2. The galactic coordi nates a nd reddening s de¬ 
termined using the m ethods of lHakkila et al.l (I1997T1 and 
ISchlegel et al.l (ll998tl 


Star 


HD 4306 
HD 26169 
HD 27928 
HD 274939 
HD 45282 
HD 108317 
HD 126587 
HD 128279 
HD 200654 
HD 218857 
BD—01° 2582 
CD—24° 1782 
CD—30° 0298 


l b 


118.057 -72.359 

289.820 -36.268 
239.682 -44.720 
254.431 -32.541 
207.050 -3.930 
286.675 67.387 
330.345 35.489 
329.070 28.356 

348.821 -41.873 
52.926 -64.414 

275.101 57.703 
217.842 -52.332 
275.051 -86.625 


E(b-v) E(b-v) 
[mag] [mag] 


0.03 

0.04 

0.02 

0.08 

0.02 

0.03 

0.03 

0.03 

0.04 

(0.82) 

0.01 

0.02 

0.04 

0.10 

0.05 

0.10 

0.02 

0.03 

0.01 

0.04 

0.00 

0.02 

0.02 

0.02 

0.01 

0.02 


ture for this subgiant would be more than 8000 K, obvi¬ 
ously an impossibly large value for a star, e.g. showing OH 
lines, which strengthens our belief that the Hakkila et al. 
method is to be preferred for our programme stars. 

The equivalent widths of the interstellar Nai doublet 
at 589 nm can also provide an independent estimate of 
the reddening; our VLT-spectra covered this wavelength 
region as well. This has been tried for a few stars for com¬ 
parison with the results of the other two methods. It is 
not always possible to make useful reddening estimates 
with this method due to blending with stellar lines or 
due to multiple interstellar components. To isolate the in¬ 
terstellar component of the lines, we performed spectrum 
synthesis using atomic dat a from the The Vienna Atomic 
Line Data Base (VALD) (iPiskunov et al.1 ll995h . As an 
example, HD 26169 shows a relatively strong main inter- 
stellar absorption a s illus trated in Figure [3 According to 
iMunari fc Zwitterl l|l997t) the equivalent width of the in¬ 
terstellar Nai D 589.0nm absorption line of 15.0pm cor¬ 
responds to E(b_ V ) = 0.05 . This is in between the red¬ 
dening estimates presented in Tabled from the Schlegel 
et al. and Hakkila et al. models. The line in the case of 
HD 45282 is much weaker (~ 2 pm) which suggests a very 
low reddening, lower than 0.05. The interstellar contribu¬ 
tion to the line will be even weaker if it turns out to be 
blended with a telluric line. Unfortunately, the interstellar 
Nai of the fast-rotator calibration star lies approximately 
at the same observer’s wavelength as the line of concern. 
Hence, it not possible to identify any telluric blend by just 
comparing the observed spectra of these two stars. 

3.3. Effective temperature: photometry 

The effective temperatures have been estimated using op¬ 
tical and IR broad-band and Stromgren photometry, in 
particular (V — K) and (b — y) colours. The V an d uvbyf3 
photometry comes from lSchuster fc Nissen ( 198 81 and pri- 
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Fig. 2. The observed (open circles) and synthetic (dashed 
line) spectrum of HD26169 around the Nai D 589.Onm 
line. Also shown is the observed spectrum of the B star 
HR 1705 (dot-dashed line), to help identifying the telluric 
lines in the observed spectrum of HD 26169. The spectral 
feature around 589.Inin is the interstellar absorption line. 
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vate communication, while the K magnitudes come from 
The Two Micron All Sky Survey (2MASS) database. The 
photometric data together with the parallaxes for the pro¬ 
gramme stars are given in Tabled As described above, we 
have used the Hakkila et al. interstellar extinction maps 
to estimate reddenings, which have been applied to de- 
redden the stellar photometry. To convert the de-reddened 
(b — y) and (V — K) colours to T„a estim ates, we have used 
the calibrations bv lAlonso et al l ( 1999 ). which are based 
on the infrared flux method: 


0 (V - K ) = 0.5558 + 0.2105(U- I\) + 0.001981(U -A') 2 - 
0.009965(U - K)[Fe/H] + 0.01325[Fe/H'] - 
0.002726[Fe/7L] 2 (1) 

0 (b _ y) = 0.5815 + 0.7263(6 — y) + 0.06856(6 — y) 2 — 

0.06832(6 - y)[Fe/H] - 0.01062 [Fe/H] - 
0.01079[Fe/H] 2 (2) 

T eff = 5040/0. (3) 

As the subgiants of the observed sample are ascending 
the red giant branch we have used the calibrations for 
giants rather than for dwarfs. We note, however, that the 
differences in temperature between the two calibrations 
are small. In our case, the dwarf calibration would have 
resulted in about 50 K lower values. 

The agreement between T e s{V — K ) and T e g(b — y) is 
quite good, as illustrated in Fig.[3] The mean difference 
is 9.70 ± 52 K (s.d.) with a maximum difference of 108 K. 
There is no obvious trend in these differences with T e g, 
although there may be a correlation with [Fe/H]. The lat¬ 
ter conclusion, however, hinges entirely on the two least 
metal-poor stars with [Fe/H] ~ —1.5. The mean of values 
























6 


A.E. Garcia Perez et al.: Oxygen abundances in metal-poor subgiants 



100 


50 

1— 

1 

0 

aT 


1 

-50 

I— 5 

-100 


4700 4800 4900 5000 5100 5200 5300 5400 
T e ,AV-K) 





100 


50 

0 

-50 


-100 


o 


- o 


o 


o 


o 

o 


o 


o O 


o 


o 


o 


O - 


-3.0 -2.5 -2.0 -1.5 

[Fe/H] 


Fig. 3. Difference in T e g when using the de-reddened 
(V — K) and (b — y) colours following Eq. and Eq. [3 
respectively. The temperature differences are shown in the 
top panel as a function of T e g (determined from (V — K)) 
and in the bottom panel versus [Fe/H]. 


based on the de-reddened (b — y )o and (V — K )o colours 
were adopted as the final stellar effective temperature. 
The resulting photometric T e g values are listed in Tabled 
together with the mean T e g without correcting for red¬ 
dening. While accounting for reddening makes a signifi¬ 
cant difference, the maximum correctio n Js nevert heless 
< 200 K when relying on the iHakkil^^^lJjllMllL^Sd- 
dening. As expected, the use of the lSchlegel et alJ 11199^1 
method would raise the de-reddened T e g values greatly in 
some cases as seen from Tab] (‘[3 As detailed below, those 
high T e ff values are in many cases inconsistent with those 
estimated independently from Ha line-profile fitting. 

The typical observational error in (6 — y) is 0.007 mag, 
which according to Eq. [2] corresponds to an error in T e g 
of about ±35 K, assuming an associated uncertainty in 
[Fe/H] of ±0.1 dex. The observational error in (V — K) 
is significantly higher, perhaps as high as ±0.05 mag, but 
due to the high sensitivity to T e g for this colour, this cor¬ 
responds to only ±45 K. These errors are consistent with 
the spread in the difference between the two colour indica¬ 
tors as shown in Fig. [3 The absolute T e g values are, how¬ 
ever, likely significantly more uncertain than this. Another 
important uncertainty stems from the adopted reddening, 
which is hardly better determined than AE(b_v) = ±0.02 
(AE(b_y) = ±0.014). This error estimate may even seem 
opt imistic in view of th e large error (0.07 mag.) quoted 
by lHakkila~ ]d] <lH>97h . However, the relatively small 


distances to our programme stars as we ll as the small 
reddening values fr om both lHakkila et al.l l|l997ll and 
ISchlegel et al.l lll998lf support such a low value. This cor¬ 
responds to an error in T e g of about ±70 K. Furthermore, 
T e ff may be affected by other systematic errors, but a rea¬ 
sonable estimate of the uncertainty in T e g is ±100 K. 


3.4. Effective temperature: hydrogen lines 


The stellar Ha 656.3 nm line is a potentially good in¬ 
dicator of effective temperature s and has been used ex ¬ 
tensively in the liter ature (e.g. iFuhrmann et al.l |l9 931: 
lHa/rklem et a,1.1 d2002ll and references therein). The fact 
that this temperature indicator is not affected by redden¬ 
ing is an advantage over photometry-based methods, es¬ 
pecially when accurate redde ning estimates are not avail¬ 
able. As discussed in detail bv lBarklem et al.l ( 2002 1 there 
are, however, several problems with this methods which 
must be overcome. An important aspect to be consid¬ 
ered in the modelling is the theory used to describe the 
broadening of the hydrogen lines, both Stark and self¬ 
broadening. Additionally, the results are dependent on 
the employed ID model atmospheres with the mixing 
length pa rameters used for com puting the convective en¬ 
ergy flux ( Fuhrmann et al . 1991 . One may worry that the 
ID model atmospheres are not sufficiently realistic to de¬ 
scribe this aspect. There is also some dependence of T e g 
as derived from H lines on the adopted stellar gravities. 
A change of 0.3 dex in logg can affect the estimates of 
the effective temperatures by as much as 60 K in the case 
of the metal-poor subgiants. Finally, the normalisation of 
the observed spectra is also a potential source of error. 
However, for our subgiants, the wings of the Ha line do 
not extend more than one spectral order so the normalisa¬ 
tion of the spectra, though not trivial is much less a critical 
issue than when dealing with turn-off stars. Furthermore, 
given the metal-poor nature of our targets, it is straight¬ 
forward to trace the hydrogen wings between the very few 
other stellar lines in this wavelength region. 


Our spectra are of sufficiently high quality in terms 
of resolving power and S/N so that the observational er¬ 
rors in the H line profiles are insignificant compared with 
the modelling uncertainties (broadening and model at¬ 
mospheres in particular). We estimate that our relative 
Ha-based T e g values are accurate to significantly bet¬ 
ter than 100 K, although the absolute values could be 
i n error by perhaps twice this amount (see discussion in 
iBarklem et al ■ 2002). As a consequence, the Ha method 
mainly serves here as an independent test of the photomet¬ 
ric values derived above, which will remain our preferred 
choice for the abundance analysis. 

The observed Ha lines are displayed in Fig.@]for five 
stars together with corresponding synthetic spectra. The 
theoretical H line profiles shown in the figure have been 
computed for the mean de-reddened T e g based on (6 — y) 
and (F — K) as well as for T e ff±100K. We note that we 
only attempt to fit the wings of the H line as the core 
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is aff ected by departures from LTE JPrzvbilla fc Butler! 
12004 and by uncertainties in the structure of the up¬ 
per atmosphere and its velocity fields. Estimates of ef¬ 
fective temperatures based on Ha are given in Tabled 
They agree with the assumed photometric stellar effective 
temperatures to within 100 K except for HD 45282 and 
CD—24° 1782. A lower effective temperature had to be 
adopted to fit their observed Ha lines, 5050 K and 4900 K 
versus 5352 K and 5228 K. The computed spectra for those 
new effective temperature values are shown in Fig.|5] 
Estimates of effective temperatures based on Ha lines 
are in this work on average about 70 K lower than the 
values based on photomet ry associated with reddenings 
from lHa.kkila et a.1.1 (4 997 1 and about 10 K higher than 
values based on the non-corrected photometry. There is no 
obvious trend in the differences of Ha temperatures and 
the photometric ones with the assumed reddening. On the 
other hand, the values of T e g coming fr om de-re ddened 
photometry using reddenings from ISchlegel et all ( 1998 1 
are higher by about 150 K than the values determined 
from Ha. This is another indication that the Schlegel et al. 
dust maps overestimate the reddenings for our programme 
stars. 


3.5. Gravities 

The stellar surface gravities were determined from 
Hipparcos parallaxes as long as the relative error A 7 t/ 7 t 
in the parallax satisfied An/n < 0.5. This was fulfilled in 
all cases except for HD 126587 and CD—30° 0298. To con¬ 
vert from parallaxes to su rface gravitie s we relied on the 
fundamental relation (e.g. INissen et a,l.lfl997 1: 


log— = log 
gO M, 


- 4 log ■ 


0 T e ff0 

0.4Vq + 0ABC + 21 og 7 r + 


(4) 


The stellar mass M was assumed to be 0.8 Mq for all 
stars which, however, is not a critical assumption for this 
me t hod. The bolometric corrections BC were taken from 
iBessell et alJ (1 9981 '). The estimates of the surface gravities 
based on Eq.QJare presented in TablcEl The main source 
of error for the stellar gravities based on this method is 
in most cases the adopted parallaxes. We estimate the to¬ 
tal uncertainty in surface gravity by adding in quadrature 
the errors due to the quoted errors in Hipparcos parallaxes 
(see Tabled, an uncertainty of 100 K in T e g and 0.05 Mq 
in mass. These ler errors are also listed in TablcEl We 
note that the effects of uncertain reddening are small as 
they affect T e g and visual magnitude in opposite direc¬ 
tions and therefore mainly cancel each other. The typical 
uncertainty in logg is ±(0.2 — 0.4). 

An independent method of deriving surface gravities 
relies on the use of theoretical evolutionary isochrones for 
the assumed mass. We have used this method as a test, 
and it is our main method for the two stars HD 126587 
and CD—30° 0298 without reliable Hipparcos parallaxes. 


We used the isochrones of iBergbusch fc VandenBerei 
A whi ch are based on the evolutionary tracks of 
IVandenBerg et TO (: 2000 b Interpolation of the isochrones 
in a-enhancement and in metallicity was performed, 
adopting 0.4 and 12 Gy for the a-enhancement and age, 
respectively, for all the stars. As corresponding theoret¬ 
ical isochrones are not available for [Fe/H] < —2.31, the 
isochrone calculated for the lowest metallicity was adopted 
for HD 126587 and CD—30° 0298 which are both quite 
metal poor, as well as for HD 200654 ([Fe/H] = —2.75). 
The rest of our programme stars have metallicities of 
about —2.3 or higher. 


The resulting surface gravities based on isochrone 
fitting are listed in Tabled Fig.|Hl shows a comparison 
of the isochrones with the surface gravities estimated 
from Hipparcos parallaxes. As can be seen, the agree¬ 
ment with the stars for which accurate parallaxes are 
available are quite good. The mean difference in log <7 
amounts to 0.21 ± 0.13. We have used these differences 
to correct the isochrone-based values for HD 126587 and 
CD—30° 0298 to place them on the same scale as the 
other stars which are based on Hipparcos parallaxes. The 
mean difference for the six stars with [Fe/H] < —2.2 is 
Alogg = 0.22 ± 0.50(s.d.). The final log <7 values for these 
two stars listed in TablcEl have thus been corrected by 
0.22. The error in logarithmic gravity for these stars is es¬ 
timated to be ±0.5. Fortunately, the method seems to be 
rather insensitive to metallicity. 


3.6. Metallicity and microturbulence 

The procedure to determine spectroscopic stellar param¬ 
eters must be iterative. The effective-temperature and 
surface-gravity estimates depend on the adopted stellar 
metallicity, which in turn depends on the other two pa¬ 
rameters. We used the Stromgren mi index to obtain a 
starting guess for the metallicity of the model atmosphere 
by which we derived new [Fe/H] values using observed 
Fell li nes. We have e mployed the same 13 Fen lines as 
in INissen et al . ( 120021 ) with the same adopted line proper¬ 
ties. These new [Fe/H] values were subsequently used to 
derive new T e g and log <7 values for the stars. New MARCS 
model atmospheres were computed with these new param¬ 
eters from which revised [Fe/H] estimates were calculated. 
The procedure was repeated until convergence which oc¬ 
curred within a few iterations. The spectroscopic metal¬ 
licity corresponding to the last iteration was adopted as 
the final stellar metallicity (Tablc^J). We note that these 
values slightly differ from those listed in TableEJ which 
correspond to the assumed model atmosphere metallicity 
for the last iteration. 

It should be noted that while knowledge about the 
stellar microturbulence is necessary in this iterative pro¬ 
cedure, the exact choice for this parameter is not impor¬ 
tant for the majority of our stars. Only for the three most 
metal-rich stars in the sample (HD 274939, HD 218857 and 
HD 45282) did the derived Fe Ii-based abundances depend 
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Fig. 4. The observed (open circles) Ha line spectra are compared for five stars with synthetic spectra (dashed lines). 
All spectra have been normalised. Three synthetic spectra are plotted for each star, corresponding to three different 
effective temperatures. The middle line corresponds to the assumed effective temperature based on the de-reddened 
photometry while the other two correspond to 100 K higher and lower T e g, respectively. The scale on the vertical axis 
goes from 0.7 to 1.1 for each panel. The observed absorption lines in the wings of the Ha line which do not have a 
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Fig. 5. The observed Ha line (circles) and the best-fitting synthetic spectra (dot-dashed lines) for HD 45282 and 
CD—24° 1782, corresponding to T e s = 5050 K and 4900 K, respectively. 


Table 4. The derived [Fe/H] values based on Fell lines 
with line-to-line scatter given. Also listed are the associ¬ 
ated errors in [Fe/H] due to uncertainties in stellar pa¬ 
rameters (100K for T e g, 0.1 dex in [Fe/H] and the errors 
in log g as given in Tabled- Finally, the impact of a mea¬ 
surement error of 0.05pm for the Fen equivalent widths 
on the derived Fe abundances is also illustrated. 


Star [Fe/H] A [Fe/H] 

AT eff Alogff A [Fe/H] A W x - A W x 


HD 4306 

-2.33 

± 

0.07 

-0.01 

0.10 

0.01 

0.07 

-0.10 

HD 26169 

-2.28 

± 

0.04 

0.00 

0.10 

0.01 

0.06 

-0.06 

HD 27928 

-2.14 

± 

0.04 

-0.01 

0.17 

0.01 

0.04 

-0.06 

HD 45282 

-1.52 

± 

0.04 

0.00 

0.05 

0.02 

0.03 

-0.02 

HD 108317 

-2.25 

± 

0.04 

0.00 

0.07 

0.00 

0.05 

-0.08 

HD 126587 

-2.87 

± 

0.05 

0.00 

0.18 

0.00 

0.07 

-0.10 

HD 128279 

-2.19 

± 

0.08 

0.00 

0.07 

0.00 

0.06 

-0.09 

HD 200654 

-2.73 

± 

0.08 

0.01 

0.12 

0.00 

0.08 

-0.11 

HD 218857 

-1.79 

± 

0.04 

-0.02 

0.14 

0.02 

0.03 

-0.03 

HD 274939 

-1.49 

± 

0.04 

-0.02 

0.13 

0.02 

0.02 

-0.02 

BD—01° 2582 

-2.11 

± 

0.07 

-0.01 

0.15 

0.01 

0.05 

-0.06 

CD-24 0 1782 

-2.25 

± 

0.08 

-0.01 

0.14 

0.00 

0.09 

-0.10 

CD—30° 0298 

-3.01 

± 

0.10 

0.01 

0.19 

0.00 

0.10 

-0.14 


on the microturbulence. The microturbulence parameters 
for those stars were estimated by requiring that the Fen 
lines should not yield any trend in Fe abundance with line 
strength. The thus estimated microturbulence parameters 


were in the range 1.1 — 1.4kms . For the other stars, the 
lines are sufficiently weak to be insensitive to the choice 
of microturbulence, which were assumed to be lffikms” 1 . 
The model atmospheres are hardly affected at all by the 
choice of the parameter. 


The equivalent widths of the Fe n lines were measured 
with the iraf task splot. Due to the weakness of the lines 
at low metallicity, all 13 Fell lines were not detected in 
all stars. In the case of CD—30° 0298 there were only 4 
lines detected while at least 11 lines were detected in the 
spectra of the rest of the stars. The choice of ionized 
iron lines was natural since these lines are less affected 
by departures from LTE (such as over-ionization) than 
neutral iron lines, even if the magnitude of the NLTE 
effects for Fel lines is still debated, mainly as a result 
of the uncertain inelastic H colli sions (compare for ex¬ 
ample the contras t ing r e sults of GraHoiiet_al . HHI; 
iThbvenin & Idiard lll999l) : iKorn et TITif^OO.'j) ) . The Fell 
lines of our prog ramme star s do not suffer from any ap¬ 
parent blends ( Nissen et IDEmi. Iron abundances for 
each of these lines were derived from th e measured equiv- 
alent widths adopting the q /-values oflNissen et al.l |2002i l 
which are based on the work bv lBiemont et al.Hl99lll . The 
mean abundance was taken as the stellar metallicity. Since 
we also performed a corresponding solar analysis using the 
same lines to derive a solar Fe abundance, the effects of 
any errors in the ^/-values on the stellar iron abundances 
relative to that of the sun are minimised. 
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Table 3. The derived stellar parameters of our sample. The effective temperatures determined from reddening 
corrected colour indexes, (6 — y) 0 and (V — K ) 0 , with reddening values according to column 4 in Tabled The means 
of those effective temperatures are presented in the fourth column denoted by T e g (adopted). In the fifth column, 
denoted by T e g (no-red), the same type of mean values are reported but with no reddening corrections applied. The 
mean values for the effective temperatures based on reddenings from ISchlegel et al.l 1 119981) are listed in the sixth 
column. In the seventh column values for effective temperature based on the Ha line profiles are listed. The remaining 
columns give values adopted for the model atmospheres for gravities (based on Hipparcos parallaxes and isochrones, 
respectively), metallicities and microturbulence parameters. 


Star 

T e ff(6 - y ) 
Hakkila et al. 

[K] 

T e s{V — K) T e ff 

Hakkila et al. adopted 
[K] [K] 

T e g T e s T e ff (Ha) 

no-red Sclilegel et al. 

[K] [K] [K] 

log g [Fe/H] £turb 

Hipparcos Isochrones 

[cgs] [cgs] [dex] [kms^ 1 ] 

HD 4306 

5021 

4958 

4990 

4890 

5024 

4940 

3.04±0.26 

2.14 

-2.35 

1.5 

HD 26169 

4962 

4982 

4972 

4906 

5181 

4972 

2.49±0.26 

2.10 

-2.29 

1.5 

HD 27928 

5040 

5047 

5044 

4976 

5078 

5044 

2.67±0.43 

2.33 

-2.16 

1.5 

HD 45282 

5390 

5314 

5352 

5208 

8170 

5050 

3.15±0.13 

3.44 

-1.53 

1.1 

HD 108317 

5301 

5300 

5300 

5263 

5338 

5300 

2.76±0.22 

3.05 

-2.27 

1.5 

HD 126587 

4694 

4730 

4712 

4590 

4907 

4812 

(1.66±0.50) 

1.44 

-2.88 

1.5 

HD 128279 

5347 

5325 

5336 

5151 

5533 

5236 

2.95±0.21 

3.16 

-2.21 

1.5 

HD 200654 

5306 

5278 

5292 

5216 

5331 

5192 

2.86F0.35 

3.01 

-2.75 

1.5 

HD 218857 

5009 

5022 

5015 

4982 

5117 

4915 

2.78±0.36 

2.39 

-1.80 

1.4 

HD 274939 

5144 

5036 

5090 

4991 

5090 

5040 

2.79±0.33 

2.73 

-1.51 

1.3 

BD—01° 2582 

5032 

5112 

5072 

5072 

5142 

5022 

2.92±0.40 

2.41 

-2.13 

1.5 

CD-24 0 1782 

5240 

5215 

5228 

5155 

5228 

4900 

3.46F0.35 

2.82 

-2.27 

1.5 

CD—30° 0298 

5175 

5216 

5196 

5158 

5234 

5044 

(2.93F0.50) 

2.71 

-3.04 

1.5 


Table 5. List of Fei lines used for iron abundances cal- 


cula tions, exc i tation al potential and log gf values from 
lO’Brian et al . (Il99lt) . 


A 

H 

X 

[eV] 

log 9f 

613.66 

2.45 

-1.41 

613.77 

2.59 

-1.35 

619.16 

2.43 

-1.42 

623.07 

2.56 

-1.28 

625.26 

2.40 

-1.77 

639.36 

2.43 

-1.58 

641.17 

3.65 

-0.72 

642.14 

2.28 

-2.01 

643.09 

2.18 

-1.95 


We have estimated the sensitivity of the derived metal¬ 
licities to the adopted stellar parameters by repeating the 
analysis with perturbed values for T e g, log g, [Fe/H] and 
£turb- For this exercise we have used our estimated uncer¬ 
tainties in the various parameters: 100 K for T e g, 0.1 dex 
in [Fe/H] and the quoted errors in logg (from Hipparcos 
parallaxes) given in Tablc|3 The changes in the [Fe/H] val¬ 
ues in each of these cases are listed in Tabled The depen¬ 
dence of temperature on the line and continuum opacity 
is such that the Fe Ii-based abundances are essentially im¬ 
mune to errors in T e ff(Alogepe < 0.02dex). Similarly, the 
resulting impact on the model atmosphere structure by 
slight changes in metallicity is very small for these metal- 
poor stars and has virtually no effect on the derived Fe 


abundances. However, the uncertainties in the stellar grav¬ 
ities lead to significant errors in the calculated iron abun¬ 
dances. In the worst cases this error may be as large as 
0.19 dex for the stars with poorly determined parallaxes, 
although in most cases the uncertainty is about 0.10 dex. 
The microturbulence influences only the iron abundance 
slightly even for the most metal-rich stars. Even a drastic 
change such as going from l.lkms -1 to 1.7kms _1 in the 
case of HD 45282 makes a difference of a mere 0.04 dex in 


[Fe/H], 

Another source of error for the Fe abundances is the 
measured equivalent widths, in particular given the weak¬ 
ness of some of the employed lines in the most metal-poor 
stars. Based on photon statistics and the measured S/N 
around the Fell lines we estimate that the uncertainties 
in the equivalent widths are at most about 0.05 pm. New 
Fe abundances with such errors added to the measured 
line strengths were derived and are also listed in Tablc^J 
Fortunately, in all cases such uncertainties in the mea¬ 
sured equivalent widths only correspond to errors in the 
derived mean Fen abundances of < 0.14dex even in the 
unrealistic case when all Fen lines are affected fully and 
in the same direction. 

In the literature and so often, one can find iron abun¬ 
dances for metal-poor solar type stars derived from neu¬ 
tral iron lines rather than from ionised lines. The neu¬ 
tral lines are easier to detect because they are in gen¬ 
eral stronger than the ionised lines. However, these lines 
may suf fer from depa r tures f rom LTE in the l i ne for - 


mation 


Grattan et all lll99fl) : iTheve nin fc Idiartl dl999l) : 


lAsulundl 12005?) : ICollet et ah 1 2005if ) so we would expect 
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Fig. 6. Comparison of the stellar gravities from Hip p arcos 
parall axes with isochrones from lBerghiisch fc VandenBerd 
1: 2001 1. The three panels display the results for three dif¬ 
ference metallicities: [Fe/H] = —2.31, —1.84 and —1.54. 
In each case the isochrone at an age of 12Gyr has been 
used. The stellar data in the [log g, logT e ff] plane are rep¬ 
resented by circles with crosses denoting the estimated 
errors. The stars are distributed between the panels ac¬ 
cording to their metallicities. 


abundances derived from them to be different from the 
iron abundances derived from Fe II lines. We investigated 
this by deriving iron abundances from the measured equiv¬ 
alent widths of a set of Fe i lines (see Table|SJ) • A trend may 
exist for abundance differences between Fe I and Fe II lines 
with the assumed metallicity; lower metallicities seem to 
show larger differences. On average, our assumed metal¬ 
licity values are 0.2 dex higher than the values derived 
from the neutral lines (see TablcEl , except for two stars 
(HD 45282 and HD 128279). Furthermore, we note that 
two of the most metal-poor stars of our sample show the 


Table 6. Stellar iron abundances determined from the 
Fei lines in TablcOO Last column reports Fel-Fen iron 
abundance differences. 


Star 

log £ Fe 

A log £ F e 
(Fei-Fen) 

HD 4306 

4.65 

-0.56 

HD 26169 

5.04 

-0.23 

HD 27928 

5.17 

-0.24 

HD 45282 

6.30 

0.26 

HD 108317 

5.24 

-0.05 

HD 126587 

4.41 

-0.26 

HD 128279 

5.39 

0.04 

HD 200654 

4.69 

-0.11 

HD 218857 

5.53 

-0.23 

HD 274939 

5.99 

-0.07 

BD—01° 2582 

5.12 

-0.32 

CD—24° 1782 

4.84 

-0.46 

CD—30° 0298 

4.22 

-0.30 


largest differences (HD4306 and CD—24° 1782). The fact 
that Fen-based metallicities are lower than the ones de¬ 
rived from Fe I lines is not easy to explain in terms of the 
NLTE effects because Fei lines should give lower values 
than Fen lines. 

Let’s now take a closer look at those objects sticking 
out of the general trend, i.e. the two stars for which the 
iron abundances derived from Fen lines are lower than 
those derived from Fe I lines, and the two metal-poor stars 
with the largest Fei—Fen differences. The first possible 
source of the observed discrepancy are the stellar parame¬ 
ters, which could have been wrongly determined for these 
four stars. In the case of HD 45282 and HD 128279, we 
note they have the highest T e g of the entire sample, as 
well as the largest assumed reddening. This may clearly 
indicate that our effective temperatures for these two stars 
have been overestimated. Indeed, a decrease of 100 K in 
their temperatures (at the border-line of our quoted uncer¬ 
tainty on this stellar parameter) would bring Fe I and Fe II 
abundances in satisfactory agreement (within 0.1 dex), 
although with opposite sign. For the latter two stars, 
HD 4306 and CD—24° 1782, we note that their Hipparcos- 
based gravities are much higher than the values derived for 
stars of similar stellar characteristics. Again, this could be 
a hint that our gravity estimates for these two stars are not 
the optimal values, which could well explain the very large 
differences observed between Fei- and Fen-based abun¬ 
dances (up to 0.5 dex!). If we were to lower their gravity 
values to, e.g., 2.5 and 3.0 respectively (from 3.04 and 3.46, 
but still within a two sigma error), we would then derive 
a Fei—Fell abundance difference of the order of 0.3dex 
which is in more accordance with the average value we 
find (0.2 dex). Note that the isochrones-based gravity val¬ 
ues for both stars are even lower than 2.5 and 3.0, which 
would lead to an even better agreement between Fe I and 
Fen. 
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W a v' e I e r g t h ( n m ) 


Fig. 7. Observed spectra with the forbidden [Oi] line 
at 630.03 nm marked with a vertical bar for four of the 
programme stars. The appearance of the spectra before 
and after correction for telluric absorption is shown; the 
corrected spectra were re-scaled along the vertical scale. 


4. Oxygen abundances 

In general it is reasonable to assume that the more differ¬ 
ent types of abundance indicators are being used, the more 
reliable are the results, provided that the various transi¬ 
tions yield consistent results. There are four different in¬ 
dicators for stellar oxygen abundances, three of which are 
located in the wavelength range covered by UVES. In es¬ 
sentially all previous work in the extensive literature on 
oxygen abundances in metal-poor stars, the abundances 
have been derived from only one or in some cases two of 
these indicators simultaneously. We have chosen to tar¬ 
get subgiants and giants near the base of the red giant 
branch at exceptionally high S/N both in the UV and 
optical regions with the aim to exploit the unique oppor¬ 
tunity to derive abundances from three basically different 
abundance indicators. 

In this section we describe how the oxygen abundances 
have been derived for the different types of lines. In the 
case of the forbidden [O i] and the OI IR triplet lines we 
have employed both the techniques of measuring equiv¬ 
alent widths and detailed spectrum synthesis, while only 
the latter is feasible for the crowded spectrum around the 
OH UV lines. The equivalent widths have been measured 
using the iraf task called splot , either by directly integrat¬ 
ing the area of the spectral line of interest or by fitting the 
spectral lines with a Gaussian or another type of function 
and then integrating that curve. As the S/N of our spec¬ 
tra is superb, the direct integration method is in general 
preferred, since a Gaussian can not perfectly describe a 
stellar spectral line due to the line asymmetries induced 
by the convective motions in the stellar atmosphere. In 
a few occasions a Gaussian fit had to be adopted due to 
presence of telluric lines disturbing the profile of the [O i] 
line. 


4.1. Solar abundances 


In order to obtain accurate [O/Fe] values for the stars, we 
must also derive consistent solar O and Fe abundances. 
The solar iron abundance was determined from each of 
the 13 Fe n lines used to derive the stellar Fe abundances 
wit h transiti on probabilities and equivalent widths from 
Nigggi^e^ , measured in the solar flux atlas of 

Kurucz et alJ lll9841 . The solar analysis proceeded com¬ 
pletely analogously to the stellar Fe determinations. The 
thus derived solar Fe abundance is log£F e = 7-56 ± 0.09 
(standard deviation). The absolute abundance is reason¬ 
a bly close to the quoted m eteoritic ab undances of 7.5 0 
in lGrevesse fc Sauvall (Il998t) and 7.45 in lAsnlundl ll2005l) . 
Indeed, had the oscillator strengths from VALD been 
adopted for these lines, solar iron abundance would have 
decreased to 7.44 while the line-to-line dispersion dimin¬ 
ished to 0.05 dex. The relatively high scatter among the 
different lines is similar t o that achieved in th e 3D anal¬ 
ysis of solar Fe II lines by lAsolund et al . (l2000t) . and pre¬ 
sumably mainly reflects errors in the oscillator strengths. 
Minor sources of error are inaccurate equivalent widths 
and departures from LTE in the line formation. Unlike 
for our metal-poor stars, in the solar case the Fen lines 
of interest are sufficiently strong (2.0-9.0 pm with 9 lines 
< 4.5 pm ) that errors in line broadening may also play 
a role. If the Unsold enhancement factor for the van der 
Waal’s broadening were changed from 1.5 to 2.0 the so¬ 
lar iron abundance would only change by 0.015 dex. Much 
more important is the impact of the microturbulence: the 
metallicity changed by —0.16 dex when the assumed value 
of microturbulence was altered from 1.15 to 2.00 km s -1 . 
By computing [Fe/H] values for each individual line for the 
stars in our sample, we removed much of the uncertainty 
in the ^/-values and could thus achieve a considerably re¬ 
duced line-to-line scatter (< 0.05 dex for seven of our stars 
and only in one case as high as 0.10dex). 

A solar O abundance was determined from the [O i] 
630.03 nm line using the m easured equivalent width of 
0.55pm from iNissen et al.l d2002ll . The 636.38nm line 
is more severely blended and was not used. As also 
the [ Oil 630.03n m line is blended, with a Nil line 
( Allende Prieto et al.l200l|l , we have subtracted its contri¬ 
bution from the f eature by adopting a solar Ni abundance 
of log evi = 6. 25 dGrevesse. fc Sauvall 199^1 and log gf = 
—2.31 llAllende Prieto et a,l.ll200ll) . The remaining equiv¬ 
alent width of 0.41pm attributable to the [Oi] 630.03nm 
line implies a solar O abundance of logeo= 8.74. That 
value was always taken as the solar O abundance indepen¬ 
dently of the type of indicator used when deriving [O/Fe] 
for the programme stars. 


4.2. Oxygen abundances from the [Oi] 630.03nm line 

Oxygen abundances were derived from the [O i] 630.03 nm 
line using both measured equivalent widths and by means 
of spectrum synthesis. The measured equivalent widths 
are given in Tablc[3 Depending on the radial velocity of 
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the star the stellar line can be blended by telluric lines 
which have to be carefully considered, in particular for 
these metal-poor stars for which the stellar line is so weak. 
As mentioned in Sect .[3 precautions were taken to avoid 
such blends and only in one case, for CD—24° 1782, the 
[O i] line turned out to be so severely blended that it was 
useless for abundance determinations. 

In six additional cases telluric lines partially blended 
the oxygen feature. In the spectra of these six stars we 
used the iraf task telluric. However, since the S/N val¬ 
ues for these B stars were lower than for the programme 
stars, this method degrades the S/N, and consequently 
should be avoided if possible. The procedure worked very 
well for HD 4306, HD 26169, HD 218857 and HD 45282, 
but less satisfactorily for HD 126587 and HD 128279. The 
residuals (i.e. the difference between the observed and the 
calculated spectra, normalised to the continuum) for these 
two latter stars were not insignificant, being of the or¬ 
der of 1-2%, so equivalent widths were not measured from 
them. Instead, a Gaussian profile was fitted to the ob¬ 
served lines to estimate the contribution of the blend to 
the wings. In the case of HD 126587 and contrary to the 
case of HD 128279, the large residuals are due to a large 
difference in airmass between the star and the correspond¬ 
ing calibration star, which makes the telluric correction 
more uncertain. Another reason was that the observa¬ 
tions were taken at different times during the night for 
the programme star and for its comparison star. It should 
be noted that in the successful cases, blending by telluric 
lines is not particularly problematic and a reliable equiv¬ 
alent width measurement could be obtained even without 
the telluric procedure. This can be seen in Fig.Q where a 
few spectra before and after applying the task are plotted. 
In the remaining stars, there was no need to perform any 
division with a B star spectrum as no telluric lines were 
disturbing the forbidden oxygen line. For those spectra, 
the equivalent widths were measured directly, using di¬ 
rect integration of the area of the feature in the observed 
spectrum. In the case of CD—30° 0298 and HD 200654 no 
convincing line detection could be made and only an up¬ 
per limit to the measured equivalent width of 0.04 pm and 
0.08 pm could be placed. The last value is based on a ten¬ 
tative measurement of the equivalent width. 

As is obvious from Tabled the measured equivalent 
widths are very small, in particular at the lowest [Fe/H]. 
It is therefore paramount to carefully estimate the un¬ 
certainties in these values. From the typical S/N of 500 
around the [O i] 630.03 nm line, we expect an observational 
error of < 0.03 pm . This calculated value is based on the 
assumption of a Poisson distribution for the signal, con¬ 
sidering noise only in the line and not in the continuum. 
Fortunately, for these metal-poor stars the placement of 
the continuum is straightforward and does not add signif¬ 
icantly to the error. This was demonstrated by perform¬ 
ing repeated independent measurements of the equivalent 
widths by different people. In all cases the rms of these 
measurements were < 0.02 pm. 


Table 7. The measured equivalent widths of the [O i] 
630.03 nm line and the resulting ID LTE oxygen abun¬ 
dances. The third column gives the oxygen abundances 
from the measured equivalent widths. The remaining 
columns describe the errors in the derived abundances 
resulting from uncertainties in T e g , log g, [Fe/H] and of 
a conservative measurement error of 0.05 pm ; the stellar 
parameter changes are the same as in Tabled] No value 
is given for CD—24° 1782 since the forbidden line is badly 
blended by a telluric line. Only an upper limit could be 
obtained for CD-30 0 0298 and HD 200654. 


Star 

BA 

[pm] 

log eo 

A log eo 

ATefi A log g A [Fe/H] AW a 

-ABA 

HD 4306 

0.25 

7.14 

0.05 

0.10 

0.02 

0.08 

-0.10 

HD 26169 

0.33 

7.05 

0.06 

0.10 

0.02 

0.07 

-0.07 

HD 27928 

0.25 

7.06 

0.05 

0.16 

0.02 

0.08 

-0.10 

HD 45282 

0.41 

7.72 

0.04 

0.05 

0.02 

0.05 

-0.06 

HD 108317 

0.25 

7.23 

0.07 

0.07 

0.01 

0.08 

-0.10 

HD 126587 

0.20 

6.32 

0.08 

0.17 

0.01 

0.10 

-0.12 

HD 128279 

0.11 

6.96 

0.07 

0.07 

0.01 

0.17 

-0.25 

HD 200654 

< 0.08 

< 6.75 






HD 218857 

0.37 

7.33 

0.04 

0.15 

0.03 

0.06 

-0.06 

HD 274939 

0.74 

7.76 

0.01 

0.11 

0.01 

0.03 

-0.03 

BD—01° 2582 

0.25 

7.17 

0.05 

0.15 

0.02 

0.08 

-0.10 


CD—24° 1782 . 

CD—30° 0298 < 0.04 < 6.38 


The [O i] line is bl ended with a Nil line at 630.04nm 
( Allende Prieto et alJ l200lfl . The theoretical equivalent 
widths for this line were computed for the stellar pa¬ 
rameters in Tabled adopting Ni abundances as the solar 
abundance logem= 6.25 scaled to the stellar metallicity. 
We h ave adopted th e qf-v alue for the Nil line given by 
lAllende Prieto et al.1 ( 2001 ). For the most metal-rich stars 
of our sample, the Ni contribution is only 0.01pm and 
much smaller for the remaining stars. Consequently, it is 
safe to assume that the entire measured equivalent width 
given in Tabled is due to oxygen. 

The measured equivalent widths have been converted 
to stellar oxygen abundances by using the LTE spec¬ 
tral line formation code eqwidth of the Uppsala spec¬ 
trum synthesis package, and MARCS theoretical 1D-LTE 
model atmospheres. The gf-v alue for the [O i] line was 
ta ken as log qf = —9.72, whi ch is the same as adop ted 
bv lAllende Prieto et al.l ( 2001 ) and lNissen et al.1 12 0021 in 
their analyses of the line. The line is in all cases sufficiently 
weak that uncertainties in the van der Waal’s pressure 
broadening, radiative broadening and microturbulence are 
completely negligible. The resulting oxygen abundances 
are given in Tabled 

We have investigated the sensitivity of the derived re¬ 
sults to the uncertainties in the stellar parameters. Tabled 
gives the oxygen abundances when perturbing the stel¬ 
lar parameters by their typical errors (±100 K in T e g , 
the quoted uncertainties in logg in Tableland ±0.1 in 
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[Fe/H]). The highest sensitivity is for log g where the er¬ 
rors are typically 0.10 dex but can be as large as 0.17 dex. 
The corresponding errors for T e g uncertainties are about 
0.05 dex while the derived oxygen abundance is insensi¬ 
tive to the employed [Fe/H]. One of the major advantages 
with using Fell lines to estimate [Fe/H] is that the sensi¬ 
tivity of the Fe II and [O i] lines to the stellar parameters 
is very similar. Indeed, the thus derived [O/Fe] values are 
practically independent of the choice of log <7 (< 0.01 dex) 
and only moderately sensitive to T e g (~ 0.05dex). It is 
therefore clear that our final error budget for [O/Fe] as 
obtained from the [ 01 ] line will be dominated by the mea¬ 
surement error. Since the abundance error will be depen¬ 
dent on AW\/W\, the estimated (conservative) observa¬ 
tional error of 0.05 pm will have the largest impact for the 
most metal-poor stars in our sample. 

Spectrum synthesis was also applied with the aim 
of checking the oxygen abundances determined from the 
measured equivalent widths. This was done using the spec¬ 
trum synthesis code BSYN which is also a part of the 
Uppsala package. The necessary line broadening due to 
macroturbulence was determined using nearby Fe 11 lines. 
The agreement between the equivalent width and spec¬ 
trum synthesis values is excellent. The equivalent width- 
based abundances were finally adopted. 

iKiselmanl lll993l) has studied the NLTE line formation 
of the forbidden oxygen line. According to his calculations, 
the LTE approximation is perfectly justified. In cool stars 
such as ours, oxygen is mostly in the ground level of the 
neutral state because of its high ionization energy and the 
high energies of the excited atomic levels. The forbidden 
line is a transition between the ground state and the first 
excited level of 01 , which are closely collisionally coupled. 
Departures from LTE are therefore insignificant for this 
line. 

4.3. Oxygen abundances from the 0 1 777.1-5 nm lines 

The Oi triplet at 777.1-5 nm falls near the edges of two 
adjacent spectral orders in the UVES standard configura¬ 
tions that we had to use due to the service mode nature of 
our observational programme. Therefore, while the S/N 
is excellent at wavelengths near the centre of the orders, 
the counts decrease rapidly towards the location of the 
triplet. Typically we measure S/N of > 250 around the 
lines but there can be substantial variations for a given 
order of a star even within the 0.5 nm wavelength range 
of the three triplet lines. An additional complication is 
residual fringing which could not be removed fully after 
flat-fielding. 

As for the [0 1 ] line, O abundances were estimated from 
equivalent widths. The oscillato r stren g ths fo r the tran¬ 
sitions were adopted from IWiese et alJ f l996 ): log gf = 
0.369, 0.223 and 0.002 for the 777.19 nm, 777.42 nm and 
777.54 nm lines, respectively. For completeness we note 
that we used the Unsold formula with an enhancement 
factor of 2.5 for the van der Waal’s pressure broadening of 



Fig. 8. Two adjacent spectral orders containing the Oi 
triplet at 777.1-5 nm, scaled to similar continuum levels. 
The positions of the triplet lines are indicated by vertical 
bars. 


the lines. This has no impact on the derived abundances 
given the weakness of the lines. Likewise, even a change 
in microturbulence by lkrns -1 introduces a change in O 
abundance from the triplet by < 0.02 dex. 

The equivalent widths of the three triplet lines were 
measured in both echelle orders for each star. In a few 
cases, no reliable equivalent width could be measured due 
to insufficient S/N or peculiar profiles. The values from 
the two orders were averaged, weighted with the S/N of 
the spectra around the lines. The quality of the spectra in 
the different spectral orders is illustrated in Fig.|SJ Based 
on the comparison of the two sets of equivalent widths, 
we estimate that the error in the mean equivalent width 
is < 0.1 pm. 

From these equivalent widths, oxygen abundances were 
derived for each line, which were subsequently averaged 
to yield the 0 1 triplet-based abundance for the stars. The 
equivalent widths and the thus derived O abundances are 
given in Tablc|S| The three lines normally give consistent 
abundances with typical line-to-line scatter of 0.08 dex. In 
one case the largest and smallest values differ by 0.22 dex. 

We have investigated the uncertainty in the derived 
oxygen abundances by varying the stellar parameters and 
the equivalent widths in the same fashion as for the forbid¬ 
den oxygen line. The results of this exercise are shown in 
Table|5| Contrary to the abundances based on the forbid¬ 
den oxygen line, those based on the triplet lines decrease 
when the effective temperatures increase due to the high 
excitation potential of the levels of these lines. The pop¬ 
ulation of the excited levels increases with temperature 
such that an increase of 100 K in T e g results in a decrease 
of 0.1 dex in abundance. As oxygen is predominantly neu¬ 
tral in the conditions prevailing in the atmospheres of our 
stars, the 0 1 lines are sensitive to the gravity. Just like for 
the [ 01 ] line, the abundance sensitivity is about 0.1 dex for 
the uncertainties of our log <7 estimates. The sensitivity to 
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Table 8. The mean and s.d. of the equivalent widths of the Oi 777.1-5 nm triplet lines measured from two different 
spectral orders and the derived oxygen abundances. The abundances given in columns 5-7 are the equivalent width- 
based LTE abundances for the three lines with the mean for those and the s.d. listed in column 8 . The last column 
give the abundances corrected for departures from LTE (see text). 


Star 


W A 




log to 




777.1 nm 

777.4 nm 

777.5 nm 

777.1 nm 

777.4 nm 

777.5 nm 


NLTE 


[pm] 

[pm] 

[pm] 






HD 4306 

0.63 ±0.16 

0.24 

0.39 

7.46 ±0.13 

7.16 

7.60 

7.40 ±0.22 

7.33 

HD 26169 

0.81 ±0.08 

0.53 ±0.10 

0.42 ±0.11 

7.39 ± 0.05 

7.32 ±0.09 

7.43 ±0.14 

7.38 ± 0.05 

7.31 

HD 27928 

0.67 ±0.06 

0.58 ±0.05 

0.28 ±0.08 

7.28 ± 0.05 

7.36 ± 0.04 

7.22 ±0.13 

7.29 ±0.07 

7.22 

HD 45282 

2.06 ±0.08 

1.77 ±0.23 

1.17 ±0.14 

7.79 ± 0.03 

7.84 ± 0.09 

7.81 ±0.07 

7.82 ± 0.03 

7.70 

HD 108317 

0.94 ±0.13 

0.69 ±0.07 

0.50 ±0.06 

7.24 ±0.07 

7.23 ±0.06 

7.29 ±0.06 

7.25 ±0.03 

7.17 

HD 126587 

0.46 

0.35 

0.15 

7.09 

7.10 

6.94 

7.04 ± 0.09 

6.97 

HD 128279 

1.07 ±0.33 

0.62 ±0.00 

0.39 ±0.11 

7.34 ±0.17 

7.21 ±0.00 

7.20 ±0.13 

7.25 ± 0.08 

7.17 

HD 200654 

0.31 ±0.05 

0.18 ±0.02 


6.74 ±0.07 

6.66 ±0.05 


6.70 ± 0.06 

6.63 

HD 218857 

1.05 ±0.13 

1.00 ±0.12 

0.57 ± 0.09 

7.60 ± 0.07 

7.71 ±0.07 

7.63 ± 0.09 

7.64 ± 0.06 

7.56 

HD 274939 

2.01 ±0.35 

1.39 ± 0.40 

0.62 ±0.23 

7.91 ±0.13 

7.81 ±0.18 

7.58 ±0.20 

7.77 ±0.17 

7.67 

BD—01° 2582 

1.03 ±0.09 

0.88 ±0.02 

0.47 ± 0.04 

7.57 ±0.05 

7.64 ± 0.02 

7.54 ±0.03 

7.58 ±0.05 

7.50 

CD—24° 1782 

0.51 ±0.10 

0.35 ±0.06 


7.27 ±0.09 

7.25 ±0.07 


7.26 ±0.01 

7.18 

CD—30° 0298 

< 0.2 



< 6.68 



< 6.68 

6.61 


log g for the triplet is fortunately largely offset by the sim¬ 
ilar effect on the Fe II lines when estimating [O/Fe]. 

The triplet is insensitive to the metallicity of the stellar 
atmosphere model. Since the triplet is quite weak for the 
relatively low T e g of our stellar sample, in particular for 
the most metal-poor stars, the derived abundances are 
sensitive to the errors in the measured equivalent widths. 
Our errors of about 0.1pm lead to an abundance error 
of as much as 0.26 dex for the weakest line in the most 
metal-poor stars, but more typically 0.05-0.10dex. 


It has since long been established that the OI triplet 
is affected by departures from LTE in late-type stars (e.g. 
iKiselmaiil l200lL and references therein). We have esti¬ 
mated the NLTE effects on the derived O abundances by 
performing specifically tailored statistical-equilibrium cal- 
cul ations fo r our st ars, in an identical fashion to those by 
iNissen et al • fc002ll . We have used a 23 level model atom 
with 43 bound-bound transition and with all 22 OI levels 
connected to O II by photo-ionization transitions. This is 
sufficiently extensive to yield reliable results since the line 
formation of the triplet is anyway largely an equivalent 
two-level problem. T he radiative da ta c omes mainly from 
the Opacity Project ( Cunto et al ■ 1993). Collisional exci¬ 
tation and ionization with electrons and charge transfer 
reactions are accounted for but inelast ic collision s with H 
have been neglected (see discussion in iKiselman l200lh ). 
The NLTE calculatio ns ha ve been performed using MULTI 
version 2.2 llCaiissonlll986i) for three different abundances 
surrounding the LTE based abundances derived from the 
triplet lines for the stars. From these calculations, NLTE 
abundance corrections were interpolated and added to the 
LTE results in TableOH For our stars, the NLTE abundance 
corrections typically amount to about —0.08 dex. 


Table 9. The sensitivity of the oxygen abundances from 
the O I triplet lines to uncertainties in the stellar parame¬ 
ters and equivalent widths. The stellar parameter changes 
are the same as in Tabled 


Star 

ATeff 

Alogc? 

A log to 
A[Fe/H] 

£ 

<1 

± 

-AWa 



LTE 



HD 4306 

-0.11 

0.11 

-0.01 

0.09 

-0.13 

HD 26169 

-0.10 

0.10 

0.00 

0.08 

-0.09 

HD 27928 

-0.10 

0.17 

-0.01 

0.09 

-0.12 

HD 45282 

-0.10 

0.04 

-0.01 

0.02 

-0.03 

HD 108317 

-0.09 

0.08 

0.00 

0.06 

-0.07 

HD 126587 

-0.11 

0.22 

0.00 

0.14 

-0.24 

HD 128279 

-0.09 

0.08 

0.00 

0.07 

-0.09 

HD 200654 

-0.09 

0.14 

0.00 

0.16 

-0.26 

HD 218857 

-0.10 

0.13 

0.00 

0.05 

-0.06 

HD 274939 

-0.11 

0.12 

0.00 

0.04 

-0.05 

BD—01° 2582 

-0.10 

0.16 

0.00 

0.06 

-0.06 

CD—24° 1782 
CD—30° 0298 

-0.10 

0.14 

-0.01 

0.10 

-0.12 


4 . 4 . Oxygen abundances from the OH UV lines 

The high quality of our observed spectra in the UV-blue 
wavelength region has also enabled the use of the OH (A- 
X) electronic lines around 310 nm for abundance determi¬ 
nation. One of the main problems with these diagnostics 
is the crowded spectral region they are located in. This ne¬ 
cessitates the use of spectral synthesis rather than equiv¬ 
alent widths for deriving oxygen abundances. The high 
S/N and resolving power of the blue spectra also allow an 
accurate placing of the continuum which is a non-trivial 
task for this region. Furthermore, it makes it easier to dis- 
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tinguish possible blends. That problem is the least for the 
most metal-poor stars in our sample where the continuum 
is relatively easy to locate and blends less disturbing, while 
the most metal-rich stars are likely to have larger abun¬ 
dance uncertainties due to this. An initial normalisation of 
the UVES pipeline reduced spectra was done by fitting the 
upper envelope of the spectra to a cubic spline function 
in IRAF. A minor re-normalisation around the OH lines 
was subsequently performed in some cases when the com¬ 
parisons between the normalised observed and synthetic 
spectra clearly required it. 

The selection of the OH lines to use was based on spec¬ 
tral synthesis for two stars bracketing the metallicity range 
of our sample, HD4306 ([Fe/H]=-2.33) and HD 274939 
([Fe/H]=—1.49). Only lines which were essentially free 
of blends and sufficiently strong also in our most metal- 
poor stars were considered. The theoretical OH wave¬ 
lengths, excitation potentials and transition probabilities 
were adopted from the line list of lOillis et ai.1 ( 2001 1. Many 
very weak lines which have no impact on the theoretical 
spectrum were culled from the original line list to save 
time in the synthesis, which was carried out with the 
bsyn code of the Uppsala synthesis package. The data 
for the atomic lines present i n this region were taken from 
VALD dPiskunov et al.ll995lf while the necessary data for 
th e few sig nificant CH lines were taken from the database 
of llvuruca ( 1993 11. The first guess for the element abun¬ 
dances used for the synthesis were the same as in the em¬ 
ployed model atmosphere, i.e. solar abundances scaled to 
the adopted metallicity with an a-enhancement of 0.4 dex. 
Individual abundances were adjusted as necessary to ob¬ 
tain the correct strengths of lines. The theoretical spectra 
were first convolved with a Gaussian appropriate to the 
instrumental resolving power after which an additional 
broadening profile, taking care of macroturbulence, rota¬ 
tion etc. was applied. The width of the latter was deter¬ 
mined individually for each star by optimising the overall 
fit between the synthetic and observed spectra both for 
atomic and molecular lines. In the end we identified seven 
spectral features between 310 and 330 nm which were suit¬ 
able for O abundance determinations in metal-poor sub¬ 
giants and giants near the base of the red giant branch. 
These are listed in TablcHIH together with data for the 
lines. Two of the features are in fact composed of two OH 
lines so data for nine lines are given. The OH 316.717 line 
is partly blended with a CH line, which, however, does not 
cause any significant problems. The O abundance deter¬ 
mined from this line is consistent with the values coming 
from the other OH lines. 

The derived oxygen abundances from each feature are 
listed in Tableful together with the mean abundance and 
standard deviation for each star. Spectral synthesis for 
all the seven lines together with the observations are dis¬ 
played in Fig.El for four stars. All lines give consistent 
results, which suggests that unidentified blends are not 
a significant problem. The typical line-to-line scatter is 
only about 0.1 dex, and 0.12 dex in the worst case. The 
observed OH lines of interest fall in general between the 


lower and upper solid lines in the figure which are the the¬ 
oretical predictions corresponding to abundance variation 
of ±0.2 dex. 

The scatter can likely be explained by the finite 
S/N , possible blends and slightly erroneous ^/-values. 
Departures from LTE could perhaps also contribute 
dHinkle & Lambert! Il975 ~) although very little work has 
been done on NLTE effects for molecular lines of elec¬ 
tronic band systems. One would expect, however, that the 
differential NLTE effects for different OH lines of similar 
excitation potential and line strength in the case of our 
lines should be marginal, even in the presence of signifi¬ 
cant overall departures from LTE in the molecular number 
densities or line-source functions. We note that there is no 
clear tendency in LTE abundance differences between re¬ 
sults from OH lines with different excitation energies or 
lines from different vibrational bands. 

We have investigated the sensitivity of the derived oxy¬ 
gen abundances to the adopted stellar parameters and dif¬ 
ferent continuum location in the observed spectra within 
reasonable limits. For the purpose, theoretical curve-of- 
growths were computed for the stellar parameters of our 
programme stars. The resulting behaviour is complex as 
the sensitivity depends both on the stellar parameters and 
the oxygen abundances. These values were perturbed with 
ATeff = 100, A log g according to the errors quoted in 
Table '[T] and A[Fe/H]= 0.1. The estimated abundance er¬ 
rors are given in Tablclwl The derived O abundance is 
quite sensitive to the adopted T e g and logg. An increase in 
T e g by 100 K would lead to about 0.2 dex higher O abun¬ 
dances while a 0.2 dex increase in logg roughly translates 
to an decrease in abundance by 0.07 dex. The sensitivity 
to the metallicity of the model atmosphere is minor. It 
should be noted that the errors due to T e g and log g goes 
in the opposite direction to those of Fen, which unfor¬ 
tunately makes the derived [O/Fe] even more uncertain. 
We have studied the sensitivity of the adopted spectrum 
normalisation in HD 45282, which has one of the highest 
[Fe/H] in our sample and thus should be the most sen¬ 
sitive to the normalisation. An alternative normalisation 
by 8%, which is on the order of what is still acceptable, 
leads to a 0.1 dex difference in the derived O abundance 
for this star. Clearly, this error is small compared with 
those associated with the stellar parameters. 

The final total errors in the OH-based [O/Fe] abun¬ 
dances were estimated for each star individually by adding 
in quadrature those from the stellar parameter uncertain¬ 
ties and the error in the mean value of the results of the 
OH feature fits. Typical resulting values are ±0.3 dex. 


5. Discussion 

5.1. Comparison between the three different indicators 

Our comparison between the results from the three oxy¬ 
gen abundance indicators is made primarily in terms of 
[O/Fe] which is a key parameter in studies of the chemical 
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Fig. 9. Comparison between observed and synthetic spectra around the lines in TablcTlTil for HD 108317, HD 126587, 
HD 128279 and HD 200654, with the synthetic spectra computed for three different abundances: the mean value for 
each star as listed in Tahl(TTTIarid these values ±0.2 dex. The position of the OH lines used are marked with vertical 
bars and their wavelengths are given in the respective panel. 

































































18 


A.E. Garcia Perez et al.: Oxygen abundances in metal-poor subgiants 


Table 11. The derived ID LTE oxygen abundances from the OH lines listed in TablelTUl The last column lists the 
derived mean abundance and its standard deviation (s.d.). 


Star 

310.327+ 

312.806+ 

312.954 

log eo 
316.717 

317.320 

318.931 

325.549 

log eo ± s.d. 


310.334 

312.810 







HD 4306 

6.65 

6.73 

6.80 

6.63 

6.68 

6.68 

6.58 

6.68 + 0.07 

HD 26169 

6.90 

7.00 

7.00 

6.90 

6.97 

6.89 

6.82 

6.93 + 0.07 

HD 27928 

6.96 

7.06 

7.08 

6.90 

6.95 

6.90 

6.86 

6.96 + 0.08 

HD 45282 

7.85 

7.90 

7.98 

7.90 

8.00 

7.90 

7.77 

7.90 ± 0.08 

HD 108317 

7.20 

7.30 

7.33 

7.25 

7.25 

7.22 

7.15 

7.24 + 0.06 

HD 126587 

6.30 

6.38 

6.40 

6.25 

6.40 

6.20 

6.22 

6.31 + 0.09 

HD 128279 

7.30 

7.40 

7.45 

7.25 

7.35 

7.25 

7.21 

7.32 + 0.09 

HD 200654 

6.82 

6.97 

7.02 

6.82 

6.82 

6.82 

6.70 

6.85 + 0.11 

HD 218857 

7.07 

7.12 

7.22 

7.07 

7.12 

7.12 

7.05 

7.11 + 0.06 

HD 274939 

7.51 

7.65 

7.70 

7.55 

7.65 

7.60 

7.56 

7.60 ± 0.07 

BD—01° 2582 

6.67 

6.82 

6.85 

6.78 

6.75 

6.75 

6.75 

6.77 + 0.06 

CD-24° 1782 

6.75 

6.95 

6.95 

6.73 

6.78 

6.73 

6.72 

6.80 + 0.10 

CD—30° 0298 

6.33 

6.41 

6.33 

6.20 

6.27 

6.16 

6.05 

6.25 + 0.12 


Table 10. Wavelengths of the selected OH lines from 
the (A-X) electronic system for oxygen abundance deter¬ 
mination, together with their corresponding vibrational 


numbers, excitation potential and log gf value. 

A 

[nm] 

v'-v" 

X 

[eV] 

log gf 

310.327 

0 0 

0.053 

-2.799 

310.334 

0 0 

0.053 

-2.355 

312.806 

1 1 

0.541 

-2.581 

312.810 

0 0 

0.210 

-2.969 

312.954 

1 1 

0.516 

-2.679 

316.717 

0 0 

1.109 

-1.544 

317.320 

1 1 

0.842 

-1.888 

318.931 

1 1 

1.032 

-1.843 

325.549 

0 0 

1.300 

-1.811 


Table 12. Errors in abundances due to uncertainties 
effective temperature of 100 K, AT e ff, in logg as given 
Table|j2 A log g , and in metallicity of 0.1 dex, A [Fe/H] 

Star 

ATefi 

A log eo 
A log g 

A [Fe/H] 

HD 4306 

0.19 

-0.06 

0.04 

HD 26169 

0.21 

-0.07 

0.03 

HD 27928 

0.20 

-0.12 

0.03 

HD 45282 

0.17 

-0.03 

0.05 

HD 108317 

0.22 

-0.07 

0.02 

HD 126587 

0.25 

-0.19 

0.02 

HD 128279 

0.20 

-0.07 

0.02 

HD 200654 

0.17 

-0.12 

0.01 

HD 218857 

0.17 

-0.07 

0.06 

HD 274939 

0.15 

-0.06 

0.06 

BD—01° 2582 

0.20 

-0.09 

0.04 

CD—24° 1782 

0.19 

-0.09 

0.04 

CD—30° 0298 

0.23 

-0.16 

0.01 


evolution of oxygen in the Galaxy. In this way, effects of 
errors in the stellar parameters are reduced in cases when 
derived oxygen and iron abundances respond in similar 
ways to changes in the parameters, i.e. for the [O i] and 
O i . 

Values of [O/Fe] and their uncertainties are sum¬ 
marised in Tabled] In Fig. m and for each indicator sep¬ 
arately, these values are shown as a function of metallicity. 
The values derived from OH (bottom panel) are by far the 
most uncertain with typical error bars in [O/Fe] of the or¬ 
der of 0.30 dex. About two times smaller than this are 
[Fe/H] and [O/Fe] error bars associated to the forbidden 
and permitted O I lines. The effective temperature plays a 
major role in the determination of accurate [O/Fe] values 
independently of the indicator used. Gravity only plays a 
significant role when oxygen abundances are derived from 
the OH UV lines. Measurement errors in equivalent widths 
also contribute to the [O/Fe] total error. 

For many decades, the accuracy of ID-LTE [O/Fe] val¬ 
ues based on [Oi] has been limited by the quality of the 
observed stellar spectra. With the era of new astronomical 
instruments, this quality has improved very considerably 
such that measurements based on the weak forbidden line, 
when not too weak, now probably give the most reliable 
estimates of [O/Fe] in metal-poor stars. Here, these esti¬ 
mates are taken as reference for the comparison between 
the indicators. 

According to Table ITdl [O/Fe] values based on OI per¬ 
mitted lines (including NLTE corrections) are on aver¬ 
age 0.19 ± 0.07 dex higher than the reference values, while 
values based on OH UV lines are 0.09 ± 0.08 dex lower. 
The magnitudes of the [O/Fe] uncertainties for the indi¬ 
vidual stars are of the same order or even higher than 
these values, except the uncertainties of the [O i]-based 
[O/Fe]. Notice that in the total error in [O/Fe], uncertain¬ 
ties associated with the analysis of the Fen lines should 
be included. These uncertainties, however, affect estimates 
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Table 13. [O/Fe] values and their uncertainties as determined from Fen lines and the three oxygen indicators: [Oi], 
OI IR and OH UV lines. The values derived from OI IR lines are based on NLTE calculations. 


Star 

[O/Fe] 

[Oi] 

AT e ff 

Alogg 

A[Fe/H] 

AHA 

[O/Fe] 

Oi 

ATeff 

Alogg 

A[Fe/H] 

AHA 

[O/Fe] 

OH 

ATeff 

Alogg 

A[Fe/H] 

HD 4306 

0.73 

0.06 

0.00 

0.01 

0.09 

0.95 

-0.10 

0.01 

-0.01 

0.10 

0.27 

0.21 

-0.16 

0.03 

HD 26169 

0.59 

0.06 

0.00 

0.01 

0.07 

0.85 

-0.10 

-0.00 

-0.01 

0.09 

0.47 

0.21 

-0.17 

0.02 

HD 27928 

0.46 

0.06 

-0.01 

0.01 

0.09 

0.62 

-0.09 

0.00 

-0.01 

0.11 

0.36 

0.21 

-0.29 

0.02 

HD 45282 

0.50 

0.04 

0.00 

0.00 

0.06 

0.48 

-0.10 

-0.01 

-0.03 

0.03 

0.68 

0.17 

-0.08 

0.03 

HD 108317 

0.74 

0.07 

0.00 

0.01 

0.09 

0.68 

-0.09 

0.01 

-0.01 

0.07 

0.75 

0.21 

-0.14 

0.02 

HD 126587 

0.45 

0.08 

-0.01 

0.01 

0.11 

1.10 

-0.11 

0.04 

0.00 

0.20 

0.44 

0.26 

-0.37 

0.02 

HD 128279 

0.41 

0.07 

0.00 

0.01 

0.22 

0.62 

-0.09 

0.01 

-0.00 

0.08 

0.77 

0.21 

-0.13 

0.02 

HD 200654 

< 0.74 





0.62 

-0.10 

0.02 

-0.01 

0.22 

0.84 

0.21 

-0.24 

0.01 

HD 218857 

0.38 

0.06 

0.01 

0.01 

0.06 

0.61 

-0.09 

-0.01 

-0.02 

0.06 

0.16 

0.19 

-0.21 

0.04 

HD 274939 

0.51 

0.03 

-0.02 

-0.01 

0.03 

0.42 

-0.09 

-0.01 

-0.02 

0.05 

0.35 

0.17 

-0.20 

0.04 

BD—01° 2582 

0.54 

0.06 

0.00 

0.01 

0.09 

0.87 

-0.09 

0.01 

-0.01 

0.06 

0.13 

0.20 

-0.24 

0.03 

CD-24 0 1782 






0.69 

-0.09 

0.01 

-0.01 

0.11 

0.30 

0.20 

-0.23 

0.04 

CD—30° 0298 

< 0.65 





< 0.88 





0.52 

0.22 

-0.34 

0.01 


based on the different criteria identically, and thus do not 
have any impact on the differences. Typical [O/Fe] error 
bars are of the order of 0.12, 0.15 and 0.3 dex for the esti¬ 
mates based on [Oi], Oi and OH UV lines, respectively. 

If our T e ff-scale turns out to be too cool, then both 
[O i] and O i-based estimates are expected to need modifi¬ 
cation. An increase of the temperatures by 100 K will bring 
the two different estimates together because the two in¬ 
dicators respond in opposite directions to these changes, 
with shifts in [O/Fe] of 0.06 dex and —0.1 dex, respec¬ 
tively. Another way of explaining the differences would 
be to advocate systematic errors in the equivalent widths. 
E.g., if we assume that the forbidden and the permitted 
lines are too strong and too weak, respectively, by their 
typical equivalent-width uncertainties as estimated above, 
the [Ol]-based [O/Fe] values should be adjusted upwards 
and the Oi-based [O/Fe] values decreased by about the 
same amount, ~ 0.1 dex, and thus reach agreement. A too 
cool T e ff-scale can also lead to underestimated [O/Fe] val¬ 
ues when based on abundance determinations from OH 
but by different amounts than for [O i]. Estimates based 
on the molecular lines increase faster with effective tem¬ 
perature than estimates based on the forbidden line (0.20 
versus 0.06 dex) and hence the two estimates get closer to 
each other with a hotter temperature scale. The mean dif¬ 
ference in [O/Fe] between values based on [O i] and OH, re¬ 
spectively, is lower than the [O/Fe] uncertainties of about 
0.20 dex due to gravity errors in the analysis of OH lines. 

Unfortunately, the standard deviations in the differ¬ 
ences are quite significant in comparison with the mean 
differences themselves. There are a few stars for which the 
differences are much higher than the mean value. Fig. llll 
even demonstrates possible relations between the differ¬ 
ences and the effective temperature; higher effective tem¬ 
perature may suggest a lower difference between abun¬ 
dance ratios from the triplet lines and the forbidden line, 
respectively. The trend for the abundance difference from 


the OH lines and the forbidden line is inverted with re¬ 
spect to this; only the coolest star, HD 126587, seems to 
break such a trend since it shows a good agreement be¬ 
tween the two indicators. Unfortunately, this is not the 
case for the triplet lines in the spectrum of the same 
star; this giant has the highest discrepancy (0.65 dex) be¬ 
tween abundances from those and from the forbidden line. 
HD 126587 has also high observational errors in O I-based 
[O/H], 0.24dex; these uncertainties together with uncer¬ 
tainties in effective temperature are, however, not high 
enough to account for the observed differences between 
results from the two different criteria. The observed equiv¬ 
alent widths of the OI IR triplet lines of this star are so low 
that ripples or fringing in the observed spectra might have 
disturbed the measurements, see Fig. El In the literature 
we can find other cases with high d iscrepancies b etwee n 
the two indicators. E.g., according to lBarbuv et all 1 2003f h 
the [O/Fe] estimates for the metal-poor giant HD 122563 
differ by 0.5 dex. 


Differences between the OH-based estimates as com¬ 
pared with those of the forbidden line are compatible with 
uncertainties in the corresponding [O/Fe] estimates if they 
are smaller than about 0.3 dex. Typically, an increase in 
effective temperature by 100 K will increase the estimates 
by 0.14 dex and reduce the differences relative to the esti¬ 
mates from the forbidden line, as will a decrease in gravity 
by an amount consistent with the expected errors in this 
quantity which will increase the estimate by about 0.2 dex. 
Also, a decrease of the equivalent width of the forbidden 
line by its maximum error, 0.05 pm, would bring down the 
[O/Fe] estimates from that line by ~ 0.1 dex. As a result 
of all this, the differences between the different [O/Fe] es¬ 
timates of about 0.3 dex would be reduced to zero. There 
are three data points, representing HD 4306, HD 128279 
and BD—01°2582 in the bottom panel in Fig.^2 which 
depart more than one would expect from the errors in 
parameters and observations. The differences plotted are 
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[Fe/H] 


Fig. 10. [O/Fe] estimates as a function of metallicity. 
Estimates based on different indicators are shown sepa¬ 
rately: [Ol] 630.03nm (top panel), Oi 777.1-5nm (middle 
panel) and OH UV 310.0nm (bottom panel). 


negative for two of them while the third one has positive 
differences (HD 128279). 

We remind the reader that our comparison between 
Fel- and Fen-based iron abundances for HD128279 may 
suggest that our derived effective temperature has been 
overestimated (cf Sect. 13. (ill . If this were the case, then the 
[O/Fe] differences plotted in Fig. El as well as the scatter 
in the [0/Fe]-vs-[Fe/H] trend would be notably reduced 
(a decrease of 100 K in T e g would lower the OH- and 
[Oi]-based [O/Fe] values by 0.17 dex and 0.04 dex respec¬ 
tively). The case of HD 4306 could be very similar, in the 
sense that for this star the Hipparcos-based gravity may be 
too large (see Sect. I3~6l> . With the newly suggested value 
of log g = 2.5, the OH-based [O/Fe] ratios would signifi¬ 
cantly increase (from 0.27 to 0.6), strongly reducing the 
abundance difference between the two oxygen indicators. 

BD—0 1° 2582 is claimed to be a CH binary star 
( Carnev et alil2003lf which may explain its departure. If 
its carbon abundance is very high then most of its oxy¬ 
gen may be in the form of CO and little in the form of 
OH in its outer atmospheric layers; thus a higher oxygen 
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Fig. 11. Departures of Ol-based (top panel) and OH- 
based (bottom panel) abundance estimates from [Oi]- 
based abundances as a function of T e g. 
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Fig. 12. UVES-spectra of HD 126587 around the Oi IR 
triplet 777.1-5 nm lines. Vertical bars point out the ex¬ 
pected positions of the three lines. 


abundance may be required to fit its observed OH lines. 
However, the [O/Fe] values based on the forbidden and 
permitted lines are also in conflict for this star. In con¬ 
trast to this star, with the worse agreement between the 
abundances from the three indicators, there is HD 108317 
with the best agreement: [O/Fe] = 0.74 ([Oi]), 0.68 (Oi) 
and 0.75 (OH). 
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We have explored whether a change of the T b q scale 
such that an assumed shift of T e g , which in turn would 
also be a function of T e g, could provide an explanation for 
the trend found for the abundance differences. We did not 
find any obvious such possibility. We noted, however, that 
the departure of our Hipparcos-based gravity estimates 
from the isochrone-based estimates depends on the value 
assumed for T e ff. The assumed gravity values of the coolest 
stars are much higher than the values based on isochrones, 
while the values for the warmest stars are much lower. 


5.2. Possible systematic errors: NLTE and 3D effects 


Up to this point, we have discussed the differences be¬ 
tween the oxygen abundances derived from the different 
abundance criteria in terms of errors in stellar parameters 
and observations. However we suspect that our 1D-LTE 
abundance estimates are affected by systematic errors. 


The analysis of the oxygen IR triplet lines was car¬ 
ried out with consideration of NLTE effects, and due to 
this the differences between the results and those of the 
[O i] line analysis were reduced. The latter i s expected to 
be immune to NLTE effects l|Kiselmanlll993 1. The [O/Fe] 
estimates may still be subject to other systematic errors, 
e.g. caused by thermal inhomogeneities in the stellar at¬ 
mospheres. At least for main-sequence stars, the agree¬ 
ment between the forbidden and permitted lines does 
not seem to improve with the use of 3D hydrodynami- 
cal mod el atmospheres , sooner the opposite seems to oc¬ 
cur llNissen et alJ2003> . Abundance corrections for oxygen 
in the 3D model a tmospheres for dwarfs of low metallic- 
ities are negative ( Nissen et SJIMi and they tend to 
be more important for [O/Fe] estimates from [Oi] than 
those from the OI triplet. There are some indi cations that 
this may be the case a lso for subgiants ( Nissen et all2002t 
(Shchukina et al.ll2004lh 3D calculations for main-sequence 
stars also suggest that the 3D [O/Fe] corrections are more 
significant at low metallicities. From Tablc^H one can see 
that the most-metal poor stars indeed tend to have the 
greatest discrepancies between the results from the differ¬ 
ent criteria. 


In main-sequence stars the OH lines are more af¬ 
fected by 3D effects than the forbidden line, leading to 
enhanced OH lines, i.e. smaller abundances. These effects 
depend mainly on [Fe/H] and to a lesser degree on effec¬ 
tive temperature; 3D corrections to [O/Fe] decrease with 
decreasin g effec ti ve temperat ure and increasing metallic- 
ity liAsplund & Oarcfa Perezll200 1i. The tendency in the 
bottom panel of Fig. HD is consistent with these computa¬ 
tions. However, we do not trace any dependence on stellar 
metallicity although 3D-effects on [O/Fe] from [O i] as well 
as from OH are both expected to depend on met allicity 
liAsplund & Garcia Perezll200l[ iNissen et al.ll2002 1. 

Schematic NLTE calculations for OH UV lines suggest 
that departures from LTE line formation occur mainly be¬ 
cause the radiation fields and sour ce functions are h igher 
than the local Planck function iAsolund & Garcia Peredl 


l2f>nih . Therefore, one would expect that NLTE abun¬ 
dance corrections for the lines in TablcHol should be 
positive (i.e. increase the abundances) and may depend 
on effective temperature as well as on metallicity. This 
may be one reason for the discrepancies between results 
based on these lines and the [Oi] line. Typical values 
f or the corrections in d warfs are of the order of 0.2 dex 
( Asolund fc Garcia PereJ200ll) and there is evidence for a 
dependency on metallicity: the corrections become smaller 
as the metallicity increases. However, these results are 
based on the two-level approach which may be not ap¬ 
propriate for OH lines, thus NLTE effects should be in¬ 
vestigated further. A full NLTE treatment is expected to 
increase the value of [O/Fe] which will compete with the 
3D-effects that will decrease [O/Fe]. 

Once 3D-hydrodynamical model atmospheres are 
available for subgiants of Pop II, 3D-effects on our [O/Fe] 
estimates should be investigated. 


5.3. The galactic evolution of oxygen - comparisons to 
other studies 


We have discussed the discrepancies between the results 
from the different indicators for the oxygen abundance 
and suggest that our (1D-LTE) results based on the [Oi] 
630.03 nrn line are to be preferred. We shall now com¬ 
ment on the galactic evolution of oxygen according to 
these results. In recent years, work done especially based 
on the forbidden line has suggested a trend in [0/Fe]-vs- 
[Fe/H] that is flat at low metallicities whi le several st ud¬ 
ies based on OH UV electronic tra nsitions lllsraelian et alJ 
j1 9981 200 it iBoesgaard et al.lll999tf suggest a linear trend, 
with [O/Fe] gradually rising to high levels at the lowest 
metallicities. 


Among our results some data points around —2.0 dis¬ 
turb the picture of a possible flat trend in [0/Fe]-vs-[Fe/H] 
at a level of 0.55 (see top panel in Fig. m- The scat¬ 
ter in [O/Fe] around that metallicity is barely significant, 
~ 0.14 (s.d.). [O/Fe] ranges from 0.41 to 0.74; the val¬ 
ues departing most from the mean correspond to HD 4306 
and HD 108317. While the last star is among those show¬ 
ing the best agreement among abundances from the three 
indicators, HD 4306 is one where the indicators disagree 
severely. HD 27928 and BD—01°2582, with low [O/Fel val¬ 
ues, are confirmed and claimed binaries ( Carnev et alJ 
booah respectively, so their abundances may be contam¬ 
inated. Note, however, tha t HD 108317 may be a binary, 
according to ICa.rnev et alJ ( 2003 1. HD 128279 has a low 
[O/Fe] value, but is another star with poor agreement 
among the indicators. On the other hand, a linear fit to 
the data points in the top panel of Fig.^Jis difficult to es¬ 
tablish, especially due to the scatter at low metallicity. We 
note that our result for the metal-poor giant HD 126587, 
with its relatively low [O/Fe] (~ 0.45), has a great effect 
on attempts to judge whether there is a ’’flat” tendency 
or a gradual increase in [O/Fe] when proceeding towards 
low metallicities. The result of a linear fit to the data 
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Fig. 13. Com parison of [O/ Fe] v a lues in Tab l dl3l with 
result s fro m: iNi ssen et a ll (1200211 . iFulbright fc Johnsonl 
1 200311) and lCavrel et alJliiToiMl . Results based on the [O i] 
line are shown in the top panel while results based on 
the O I triplet are shown in the bottom panel. A consis¬ 
tent value for the solar oxygen abundance, 8.74 dex, was 
adopted for all the data points. 


points is [O/Fe] = -0.09(±0.08)[Fe/H] + 0.36(±0.15). If 
instead, metallicities based on the Fei are assumed, this 
trend becomes steeper. This is due to the dependence on 
metallicity of the iron abundance differences discussed in 
Sect. 13.61 


We have chosen to compare our results with 
four studie s in which th e met all icity ranges overla p 


with ours: llsraelian et al. 1 

1998|); Nissen et al. (2002|); 

Fulbright & Johnson (20oj); 

Cavrel et al. (20041). The re- 

suits in 

INissen et al. (|2002h 

are easy to compare with 


because we have used the same method and th e same 
tool s. The m etho ds use d by IFulbright fe Johnsonl (120031) 
and ICavrel et alJ (120041) are also very similar to ours al- 
thoug h the tools are par tly different. The method utilised 
bv llsraelian et al ., ( 1199 8?) is somewhat deviating from ours. 
No attempt was made to homogenize the scales of stel- 
l ar parameters in the com parison with them. In the 
IFulbright fc Johnsonl d2003i) case, we could use the option 
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Fig. 14. Compariso n between our results and those of 
llsraelian et al.l (1998) as regards [O/Fe] values based on 
the OH line at 316.7 nm. The bottom panel shows how 
the [O/Fe] values based on the single OH line at 316.7 nm 
deviate from the mean values for all OH lines used in the 
study. The same solar oxygen abundance, 8.74, was used 
for all data points. 
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to compare to abundance results based on the Alonso et 
al. temperature scale. [O/Fe] values plotted in Fig.Udlwere 
calculated from the values of [Fe/H] and log eo as reported 
in the different papers after corrections of oxygen abun¬ 
dances for differences in the assumed log (//-values. For the 
solar oxygen abundance, we took 8.74 dex, see Sect. 14. II 

The different studies cover d if ferent stellar evo- 
lutionary p h ases. While llsraelian et al.l (Il998tl and 
Nissei^etjd (1200211 mainly stud y main-sequence stars, 
Fulbright, fc Johnsonl (120031) and lCavrel et al.l (l2004h anal¬ 
yse subgiants or giants. Our study goes to lower metallic¬ 
ities than the studies of Nissen et al. and Fulbright et 
al. but not as low as Cayrel et al. (Fig. El- In general, 
our [O/Fe] values, whether based on the forbidden or the 
triplet lines, follow the same tendencies with [Fe/H] as in 
the other studies. For oxygen abundances from the for¬ 
bidden line, see top panel in Fig. El there is no clear dis¬ 
tinction between our and the rest of the data points. For 
the triplet line results (bottom panel) there is a system¬ 
atic difference between, on one hand our results and those 
of Nissen et al., relative to those of Fulbright et al. on 
the other; a differentiation which seems to occur between 
dwarfs and giants. In that plot, most of our subgiants lie 
close to the dwarfs from Nissen et al., especially the not 
very cool ones. A few of the cooler stars get instead close 
in [O/Fe] to the giants of Fulbright et al. 

We have three sta rs in common with 
IFulbright fc .Johnsonl <20031) : HD 45282, HD 108317 
and HD 218857. These authors derived oxygen abun¬ 
dances from both the permitted and the forbidden oxygen 
lines using two different stellar parameter scales. We 
have compared our results with their results for the 
Alonso et al. scale with consideration of the differences 
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in the assumed stellar parameters, and find that our 
and their abundances agree to within 0.1 dex. We note 
that they use Kurucz model atmospheres and that their 
oxygen abundance estimates are based on the use of both 
forbidden lines. 

In order to com pare our results from O H with the lin¬ 
ear fit suggested bv llsraelian et al. (Il99 8f). it is necessary 
to displace their fit by 0.2 dex along the vertical axis to ac¬ 
count for the difference in the adopted solar oxygen abun¬ 
dance. Our results lie below theirs. This could be due to 
the use of a different set of OH UV lines or due to the use 
of different oscillator strength values. These possibilities 
have been explored by a comparison based on the only OH 
line in common (316.7 nm). The oxygen abundances were 
corrected for the differences in the logg/-value (—1.694 
versus —1.544), and plotted in Fig.^] Obviously, our data 
points lie significantly below those of Israelian et al.; only 
the subgiants with the highest T e g may be close. We re¬ 
mind the reader that the effective temperature for two 
of these stars (HD 45282 and HD 128279) may have been 
overestimated and that if a decrease of 100 K in their tem¬ 
peratures (see Sect. rr?ii is taken into account the data 
points corresponding to these stars will move down in the 
figure and consequently away from the ones of Israelian 
et al. However, the decrease in gravity suggested for other 
two stars of our sample (HD 4306 and CD—24° 1782) in 
Sect. mH would instead move up these two objects. In con¬ 
clusion, although star-to-star differences may change if the 
above mentioned changes in T e s and log g respectively are 
applied, the mean difference between our abundances and 
Israelian et al. is almost unaffected. 

The bottom panel in Fig.ll4lshows how the [O/Fe] val¬ 
ues based on the single line depart from the mean val¬ 
ues for all OH lines used in each study. These depar¬ 
tures are quite small compared with the differences be¬ 
tween us and Israelian et al. Those differences are prob¬ 
ably due to the fact that Israelian et al. analysed stars 
close to the turn-off point. These stars are expected to 
show greater effects due to thermal inhomogeneities on 
the ab undances derived from OH lines than our cooler sub - 
giants JAsplund fc Garcia Perezl200ltlNissen et al.l2002if . 

In general, the systematic differences in oxygen abun¬ 
dances between dwarfs and giants found in the present 
study may most probably be explained as results of a 
change with effective temperature and surface gravity of 
the systematic errors, as was already suggested by line 
comp utations for O H UV lines in 3D model atmospheres 
by Asplund fc Garcia Perez) ( 20011) . 


6. Conclusions 

Until 3D hydrodynamic model atmospheres for subgiants 
are available which will make 3D-NLTE line formation 
computations possible, 1D-LTE analyses based on the 
forbidden 630.03 nm line of oxygen will probably pro¬ 
vide the best estimates of stellar [O/Fe] when observ¬ 
able. These values are by far the most accurate, with typ¬ 
ical errors of the order of 0.10 dex. Note although that 


they may not be free of significant systematic errors as¬ 
sociated with inhomogeneities in the stellar atmospheres. 
Our best [O/Fe] estimates suggest a value at a level of 
0.55 ± 0.13 for Pop II subgiants, relatively independent 
on metallicity; the best linear fit to these estimates is 
[O/Fe] = — 0.09(±0.08)[Fe/H] + 0.36(±0.15). 

Both [O/Fe] values based on the IR triplet lines and 
on the OH UV lines depart from the the [O i]-based esti¬ 
mates. The departures from the [O i] results are more sig¬ 
nificant for the abundances derived from the triplet lines, 
on average 0.19 ± 0.22 dex (s.d.), than for the values from 
the molecular lines, —0.09 ± 0.25dex (s.d.). In the case 
of the triplet lines, the departures decrease with increas¬ 
ing effective temperatures while in the case of the OH 
lines they increase, turning from being negative to posi¬ 
tive. According to line computations for 3D model atmo¬ 
spheres of main-sequence stars, 3D abundance corrections 
are negative and decrease with decreasing effective tem¬ 
perature which is not in contradiction with the present 
observations. In general, 3D effects on abundances and de¬ 
partures from LTE depend on stellar parameters and this 
can be the reason for the fact that differences in [O/Fe] 
between the indicators depend on the fundamental stellar 
parameters. 

In order to increase the accuracies of [O/Fe] the S/N 
in the observations should be further improved and the 
effective temperatures should be determined more accu¬ 
rately. Accurate gravities are as well necessary to improve 
the values based on OH UV lines whose typical errors are 
of the order of 0.3 dex. 

Our results for subgiants follow the trends found by 
many others for main-sequence stars and giants in the 
last few years with the exception of our OH-based results 
which are significantly lower. 
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